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Physical and chemical systems can be characterized by 
their natural frequency uj and energy scales. It is hardly 
an exaggeration that most of what we know about such 
systems, from the acoustics of a violin to the energy lev- 
els of atoms, comes from their response to perturbations 
at these natural frequencies. For instance, chemists and 
biologists use infrared spectroscopy absorptions around 
1650 cm^^ (which corresponds to frequencies of 49 THz, 
wavelengths of 6 /xm, and energies of 0.2 eV) to identify 
the carbon-carbon double bond in organic compounds. 
And it was the observation of light emission from atoms 
at discrete energies in the electron Volt (eV) range that 
led to the quantum theory. 

It is of course the same situation in 'correlated' elec- 
tron materials. We can learn about the novel effects of 
strong electron-electron interactions and the properties 
of collective states of matter (superconductors, quantum 
magnets etc.) by characterizing their response to small 
amplitude perturbations at their natural frequencies. In 
solids, these natural frequency scales span an impres- 
sively large frequency range from x-ray down to DC. This 



incredibly broad range means that a blizzard of experi- 
mental techniques and analysis methods are required for 
the characterization of correlated systems with optical 
techniques. 

This short review and lecture notes attempt to lay out 
a brief summary of the formalism, techniques, and anal- 
ysis used for 'optical' spectroscopies of correlated elec- 
tron systems. They are idiosyncratic, occasionally opin- 
ionated, and - considering the breadth of the subject - 
incredibly brief. Unfortunately, there is no single com- 
plete treatise that presents a complete background for 
these topics in the context of correlated electron materi- 
als. However, there are a number of excellent resources 
that collectively give a solid background to this field. 

I recommend: 

1. F. Wooten, "Optical Properties of Solids", (Aca- 
demic Press, New York, 1972). 

The truly excellent and classic introduction to the 
subject of the electrodynamic response of solids (al- 
though there is very little in it in the way of elec- 
tronic correlations). The book that multiple genera- 
tions of spectroscopists learned from either directly 
or indirectly. To say that it is full of only basic 
obvious things is akin to complaining that Mac- 
beth is full of quotations. Unfortunately it is out 
of print. . . . ([I]) 

2. M. Dressel and G. Griiner, "Electrodynamics of 
Solids: Optical Properties of Electrons in Matter" , 
(Cambridge University Press, 2002). 

Much newer and modern with many plots of rele- 
vant response functions. Good discussion of many 
important techniques HP. 

3. G. D. Mahan, "Many-Particle Physics", (Plenum, 
2nd ed., 1990). 

The go-to resource for pertubative treatments of 
correlations in solids (3^. Excellent discussion of 
sum rules. 

4. Richard D. Mattuck, "A Guide to Feynman Dia- 
grams in the Many-Body Problem", (Dover, 2nd 
ed., 1992). 

It has been called "Feynman diagrams for dum- 
mies." Well ... I like it anyways (4J- 

5. M. Tinkham, Introduction to Superconductivity, 
2nd Ed. , McGraw-Hill, NY, (1996). 

The classic exposition on superconductivity. Ex- 
cellent treatment of EM response of superconduc- 
tors (and other BCS-like states) both from the phe- 
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nomenological two-fluid model and the full Matthis- 
Bardeen formalism (5). 

6. E. van Heumen and D. van der Marel, "Salerno 
lectures: Optical probes of electron correlations 
in solids" , Lecture notes for "XI Training Course 
in the Physics of Strongly Correlated Systems", 
Vietri sul Mare, Italy, 2-13 October 2006. Pub- 
lished in Lectures On The Physics Of Strongly 
Correlated Systems XI, Editors F. Mancini and 
A. Avella; Available at http : //optics . unige . ch/ 
vdm/marel_f iles/ salerno .lectures .pdf 

Very good lecture notes some of which parallels the 
treatment here (6). 

7. D. van der Marel, "Optical signatures of electron 
correlations in the cuprates" , "Strong interactions 
in low dimensions. Series: Physics and Chem- 
istry of Materials with Low-Dimensional Struc- 
tures" Vol. 25, 237-276 (2004); Also available at 
[http : / / arxiv . org/ abs/ cond-mat/0301506 
Excellent treatment of underlying theory, particu- 
larly that associated with various sum rules. Many 
relevant examples 

8. Andrew Millis, "Optical conductivity and corre- 
lated electron physics" , "Strong interactions in low 
dimensions. Series: Physics and Chemistry of Ma- 
terials with Low-Dimensional Structures" Vol. 25, 
195-235 (2004); Also available at http://phys. 
Columbia. edu/~iiiillis/july20 .pdf 
Excellent review covering many aspects of the the- 
ory of optical measurements for correlated electron 
physics. Particularly good discussion on the use of 
sum rules. Explicit calculations within Kubo for- 
malism. Many experimental examples with accom- 
panying theoretical discussion. Unique treatment 
on the applicability of the conventional descriptions 
of the coupling of EM radiation to solids ^ 

9. D. Basov and T. Timusk, "Electrodynamics of 
high-Tj, superconductors". Reviews of Modern 
Physics 77, 721 (2005). 

A thorough review of the use of optical probes in 
the cuprate superconductors (0). 

10. L. Degiorgi, "The electrodynamic response of 
heavy-electron compounds" , Rev. Mod. Phys. 71, 
687 (1999). 

11. A. Millis and P. A. Lee, "Large-orbital-expansion 
for the lattice Anderson model", Phys. Rev. B 35, 
3394 (1987). 

Experimental and theoretical papers with excellent 
sections on the phenomenological expectations for 
the optical response of heavy fermion materials and 
optical self-energies 

12. C.C. Homes, "Fourier Transform Infrared 
Spectroscopy", Lecture notes available here 
[http:/ /infrared. phy.bnl.gov/pdf/homes/fir. pdf 

A good introduction to various technical issues as- 
sociated with Fourier Transform Infrared Reflectiv- 
ity (FTIR); the most commonly used measurement 
technique for optical spectra of correlated electron 



systems JiWjl . 

13. George B. Arfken and Hans J. Weber, "Mathe- 
matical Methods for Physicists" (Academic Press, 
2005). 

The classic textbook of mathematical methods for 
physicists. Now with Weber. Good reference for 
Kramers- Kronig and Hilbert transforms (13). 

14. John David Jackson, "Classical Electrodynamics" 
(3rd ed., Wiley, 1998). 

No motivation need be given 

Other references as cited below. 



I. INTRODUCTION 

As mentioned above, the energy and frequency scales 
relevant for correlated systems span the impressive range 
of the x-ray down to DC (Fig. [ij. For instance, atomic 
energy scales of 0.5 eV to the keV make solids possi- 
ble through chemical bonding. These energies manifest 
themselves explicitly in correlations by, for instance, set- 
ting the energy scale of the large on-site repulsions C/ (1 
- 10 eV) of electrons which can lead to Mott-Hubbard 
interactions and insulating states. Typical overlap inte- 
grals between atomic wavefunctions in solids are at the 
low end of the eV energy scale and set the scale of Fermi 
energies in metals and hence the energy scale for electron 
delocalization and roughly that also of plasmon collective 
modes. Various other collective modes such as phonons 
(lattice vibrations) and magnons (magnetic vibrations) 
are found at lower energy scales, typically at fractions 
of an eV. At energies of order 50 meV are the super- 
conducting gaps of optimally doped cuprates. At even 
lower energies of the few meV scale the scattering rate of 
charges in clean metals is found. This is also the energy 
scales of gaps in conventional superconductors. Many lo- 
cal /-electron orbitals, which are relevant in Kondo ma- 
terials are found at these energies. At even lower scales 
(approximately lO's of /leV) can be found the width of 
the 'Drude' peaks found in the AC conductivity and the 
superconducting gaps of clean heavy fermion systems. 
Finally the 1-10 /ieV energy scale corresponds to the 
lowest temperatures typically reached with conventional 
cryogenics in the solid state (10 - 100 mK). The corre- 
sponding frequencies become important when studying 
crossover phenomena near quantum critical points with 
hu) K, ksT. All such energy scales can be studied in a 
number of different contexts with various photon spec- 
troscopies. 

Although there are many different photon spectro- 
scopies that can be discussed, that span these scales, in 
these lectures I concentrate on 'optical' spectroscopies, 
which I define as techniques which involve transitions 
with net momentum transfer q = and whose absorp- 
tion and polarization properties are governed by the op- 
tical dipole matrix element. I do not discuss the fas- 
cinating and important work being done using other 
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FIG. 1 (Color) The electromagnetic spectrum from radiofrequencies to gamma rays. Diagram from the LBL Advanced Light 
Source web site http://www.lbl.gov/MicroWorlds/ALSTool/EMSpec/EHSpec2.html. 



photon spectroscopies using light and charge in corre- 
lated systems with for instance, Raman spectroscopy 
(flSjl . electron energy loss, Brillouin scattering, optical 
Kerr rotation, photoeniission, and fluorescence spectro- 
scopies to give a tremendously very incomplete list. In 
the present case, the quantities of interest are typically 
the complex frequency dependent conductivity a{uj) or 
dielectric constant e(ti'). Alternatively data may be ex- 
pressed in terms of one of a variety of straightforward 
complex parametrizations of these quantities like the in- 
dex of refraction n and absorption coefficient k, or the 
complex surface impedace Zg. These quantities are de- 
fined for the interaction of light with materials generi- 
cally from zero frequency to arbitrarily high frequencies, 
in these lecture notes, however, I confine myself to the 
frequency range from microwaves, through the THz, to 
infrared, visible and ultraviolet. Each of these regimes 
gives different kinds of information and requires differ- 
ent techniques. For instance, microwaves are measured 
in cavities, striplines, or with Corbino techniques. THz 
is a huge growth area with the advent of time-domain 
THz and the increased use of Backward Wave Oscillators 
(BWOs) . The infrared and visible range are measured by 
Fourier Transform Infrared Reflectivity (FTIR). Visible 
and Ultraviolet can be measured by grating spectrome- 
ters and spectroscopic ellipsometry. Of course there are 
large overlaps between all these regimes and techniques. 

I caution that throughout these lecture notes, I fre- 
quently make use of the language of quasi-free electrons 
with well-deflned masses and scattering times. It is not 
clear that such a description should be automatically 
valid in correlated systems. In fact, it is not even clear 



that such a description should be valid in 'normal' ma- 
terials. Why should an ensemble of lO^'^ electrons/cm^ 
interacting with each other via long range Coulomb in- 
teraction have excitations and states that resemble any- 
thing like those of free electrons? Naively one would ex- 
pect that excitations would be manifestly many-electron 
composite objects wholly unrelated to the individual par- 
ticles which constitute the system. It is nothing short of 
a miracle that, in fact, many materials are relatively well 
described by the assumption that interactions do not play 
a principle role in the explicit physics. Low-energy ex- 
periments (e.g. DC resistivity,heat capacity) on materials 
like sodium, gold, or silicon indicate that many aspects 
of such systems can be well described by free-electron 
physics. The largest affect on electrons appears to be the 
static fleld of the ionic cores and the average static effect 
of the other electrons. In other materials, like for instance 
the heavy fermion compounds, experiments seem to in- 
dicate that the charge carriers are free, but have masses 
many times larger than that of a bare electron. But even 
here at the lowest energy scales there are no explicit signs 
that charger carriers are simultaneously interacting with 
10^'^ other charges and spins. In general, it is surprising 
that considering the close proximity that many electrons 
have to each other that such interactions appear to be 
only quantitative, but not qualitative perturbations on 
the free-electron physics. 

In a system of interacting fermions, the relatively weak 
effect that even potentially strong interactions have on 
the underlying physics can be understood by realizing 
the strong constraint that Fermi statistics and a Fermi 
surface provides on the scattering kinematics. This phase 
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space constriction gives, within the conventional treat- 
ment, a scattering rate that goes Hke + at low 
frequency and temperature. Due to this quadratic depen- 
dence, for some small frequency the particle's scattering 
rate will be less than its energy and one can say that the 
quasi-particle excitation is well-defined and scattering is 
a minor perturbation on the free-electron physics. Such 
a perturbation renormalizes to zero in the w — > limit. 
The effect of interactions can be subsumed into giving 
quasiparticle excitations a finite lifetime and a renormal- 
ized energy parameterized by an effective mass m*. 

Landau hypothesized that if one envisions slowly turn- 
ing on an interaction potential in a gas of non- interacting 
electrons, some character of the original system would re- 
main ([ro]). Specifically, he conjectured that there would 
be a one-to-one correspondence between states and ex- 
citations of the non-interacting system and those of the 
interacting system. In which case there can be said to be 
an adiabatic continuity between the two and one can try 
to understand and model the interacting system by mod- 
eling the non-interacting one. A system with such a map- 
ping is termed a Fermi-liquid. Its quasi-free electron- 
like excitations are termed quasi — particles. 

The success of Fermi liquid theory comes from, as 
mentioned above, the constraint on scattering kinemat- 
ics for low-energy excitations. This gives the result that 
quasi-particle excitations are only defined at arbitrar- 
ily small energy scales. The condition for well defined 
quasiparticles to exists, <C w, means that as one in- 
creases the parameter /3, which characterizes the strength 
of electron-electron excitations, the maximum energy of 
well-defined excitations decreases, but well-defined quasi- 
particles will still exist at low energies. Within this view, 
stronger interactions mean only that Fermi liquid behav- 
ior will occur at lower-energy scales i.e. at lower temper- 
atures and excitation energies. 

Systems for which the Fermi-liquid paradigm is valid 
can be described at low energies and temperatures in 
terms of quasi-free electrons. Although the vast ma- 
jority of normal metallic systems do seem to obey the 
Fermi liquid phenomenology, it is unclear whether such 
a description is valid for 'correlated' electron systems. 
A violation could occur for instance in Hubbard models 
when the intersite hopping parameter t, which parame- 
terizes the band width, is much smaller than the onsite 
Coulomb interaction U resulting in an insulating state. 
So far however the fractional quantum Hall effect is the 
only case where Landau's conjecture regarding a one-to- 
one correspondence of states has been experimentally 
falsified (jTTI IT^ . On the theoretical side, exact solutions 
to ID interacting models also show definitively that the 
fundamental excitations in ID are not electron-like at all, 
but are fractions of electrons: spinous and holons that 
carry spin and charge separately !TT)| . This would be an- 
other case where the electronic quasiparticle concept is 
not valid, but thus far in real systems the residual higher 
dimensionality has been found to stabilize the Fermi liq- 



uid (See for instance Ref. (I5U|) ^ 

It is true that many systems that we call strongly cor- 
related exhibit a phenomenology inconsistent with the 
Fermi-liquid paradigm. It is unclear whether this is be- 
cause interactions have driven such materials truly into 
a non-Fermi liquid state or because the Fermi liquid phe- 
nomena is obscured at the experimentally accessible tem- 
peratures and frequencies. For instance, it is important 
to keep in mind that the much heralded non-Fermi liquid 
behavior of high Tc cuprates is actually exhibited at rel- 
atively high temperatures (^^100 K) above the occurence 
of superconductivity. It may be that the inopportune oc- 
currence of superconductivity obstructs the view of what 
would be the low-energy quasiparticle behavior. 

At this point it is far from clear to what extent the 
Fermi liquid paradigm is valid in many of the materials 
being given as examples below. In this regard, one must 
keep in mind however, that despite the murky situation 
much of the language tossed around in this field presup- 
poses the validity of it and the existence of well-defined 
electronic excitations. Indeed, much of the used terminol- 
ogy shows this bias. The terms 'density of states', 'effec- 
tive mass', 'scattering rate', 'Pauli susceptibility', 'band 
structure', 'electron-phonon' coupling all require the con- 
text of the Fermi liquid to even make sense. Clearly such 
language is inappropriate if such excitations and states 
do not exist! Despite this, in the literature one sees many 
papers discussing, for instance, optical or photoemission 
electronic self-energies, while at the same time the au- 
thors discuss the non-Fermi liquid aspect of these mate- 
rials. It is not clear whether it is appropriate to discuss 
electronic self-energies in materials where the elementary 
excitations are not electron- like. It is the case that while 
the formalism for generating electronic self-energies from 
optical or photoemission data as discussed below may 
be followed straightforwardly, the physical signficance of 
such self-energies is not necessarily clear. For instance, 
although the parametrizations of optical spectra in terms 
of the frequency dependent mass and scattering rate from 
the extended Drude model (see below) can always be 
valid as a parameterization, it is reasonable that one 
can only assign physical significance to these quantities 
if the quasi-particle concept is valid in the energy range 
of interest. 

Some aspects of the below formalism arc model in- 
dependent and some rests on the concept of well-defined 
electronic Fermi-liquid excitations. Although I will try to 
make the various distinctions clear, I use the language of 



^ Carbon nanotubes, which are in principle perfect ID systems, 
have been claimed to have tunneling conductances that obey a 
voltage - temperature scaling that is consistent with Luttinger 
(ID) liquid physics (i21l) . However there are serious questions 
whether this study showed scaling to the requisite precision in 
order to confirm a non-Fermi liquid nature. It may also be pos- 
sible to account for the data in terms of temperature dependent 
Coulomb blockade (l22ll . 
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quasi-free electrons almost entirely throughout the below, 
because at the very least it provides a rough intuition of 
the kinds of effects one expects in insulator, metals, and 
superconductors. It also provides a self-contained formal- 
ism for the analysis of optical spectra. A generalization 
of these ideas to strongly correlated systems does not cur- 
rently exist. In the spirit of learning to walk before one 
learns to run, I use the language of quasi-free electrons 
in these lecture notes throughout, but I caution on the 
naive application of these ideas, which are only formally 
true for non-interactings systems to strongly interacting 
ones! The generalization to non-Fermi liquid, correlated, 
strongly interacting etc. etc. systems is left as an excer- 
cise for the reader!!! 



II. FORMALISM 

A. Optical Constants of Solids 

Any discussion of the interaction of light with matter 
starts with Maxwell's equations 



P = XeE, 
M ^ XrnH, 

J = aE. (9) 

Typically we express the electric and magnetic suscep- 
tibilities in terms of dielectric functions e = 1 -I- 47rxe and 
magnetic permittivity = 1 + 47rx,„. Except for explic- 
itly magnetic materials, Xm = 0. The dielectric function 
(or its Fourier transform) is a response function that con- 
nects the field E at some time t and position r with the 
field D at some later time and position. Generally we 
first define it in the time and position domain via the 
relation 

D{r,t)^ I I t{r,r',t,t')E{r',t')Sr'dt'. (10) 

One could also describe the system's response in terms 
of current and conductivity (t(w) 



V ■E{r,t) = A^p{r,t), (1) 

V xE{r,t)^-^-^B{r,t), (2) 
c ot 

^■Bir,t) = 0, (3) 

WxB{r,t)^-lE{r,t) + -J{r,t), (4) 
c ot c 

where E and B represent the electric and magnetic fields 
averaged over some suitable microscopic length (typically 
the incident light wavelength). As usual we introduce 
auxiliary fields as 



D = E + AttP, 
H = B- 47rAf, 



(5) 
(6) 



where D and H have their conventional definitions and 
M and P are magnetization and polarization. Using the 
continuity equation for the electric current V • J = — ^ 
and accounting for all sources for free conduction, po- 
larization, and magnetization currents Jtot = Jcond + 
^ -I- cV X Af as well as the external and induced charges 
Ptotai = Pext+Pind we get the additional equations (please 
refer to Jackson (fT4|l for a more extended discussion). 



J(r,i) 



a{r,r\t,t')E{r',t')(fr'dt'. (11) 



In this context, the real space and time dependent e 
and a are usually referred to as memory functions for ob- 
vious reasons. Note that the principle of causality is im- 
plied here; effects cannot precede their causes. For anal- 
ysis of optical spectra we are typically more interested 
in their Fourier transforms a^q^ui) and e{q^uj). This is 
because these are the quantities which are actually mea- 
sured, but also because the mathematical complications 
of convolution in Eqs. [TO] and [iT] is avoided; convolution 
in the timed-domain becomes multiplication in the fre- 
quency domain. The conductivity and dielectric function 
are related by^ 



(12) 



Given the very general form of Eqs. [TO]and[TT| simple 
physical considerations allow a number of general state- 
ments to be made. First, due to the vast mismatch be- 
tween the velocity of light and the typical velocity of 
electrons in solids, we are typically concerned with the 
q = limit of their Fourier transforms ^ . This means that 
except in a few circumstances, where one must take into 



IdD An 



J P.: 



in 



-J. 



cond ■ 



(7) 

(8) 



Here we assume that we are in the linear regime and 
the response of polarization, or magnetization of current 
is linear in the applied field. We therefore write 



One frequently sees the complex response function written using 
the real part of the conductivity and the real part of the dielectric 
function e.g. e = ei -|- lAirai/uj. 

Alternatively one may say it is the q = limit is relevant be- 
cause of the large difference between the momentum carried by 
an 'optical' phonon and typical lattice momenta 
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account non-local electrodynamics, (very pure metals or 
clean superconductors for example), one can assume that 
down to the scale of some microscopic length there is a 
local relationship between the quantities given in Eqs. |9] 
This means that while these quantities may have spa- 
tial dependence (for instance the current J is confined to 
surfaces in metals), the proportionality expressed in Eqs. 
[9] holds. One can use the above expressions to rewrite 
Maxwell's equations explicitly in terms of e, ^ and a. 
(Please see Jackson (fT4| for further details.) 

Furthermore, the principle of causality - that effects 
can not proceed their causes - demands strict temporal 



considerations regarding the integrals in Eqs. 10 and 11 
This leads to the powerful Kramers-Kronig relations re 
lating the real and imaginary part of such response func 
tions. We can rewrite Eq. [TT] as 



J{t) 



a{t~t')E{t')dt'. 



(13) 



From causality the conductivity memory function a(t- 
t') has the property that a{T < 0) = with t ^ t — t' 
The Fourier transform of a{t — t') is then 



(t{u)) 



CT{T)e''^^dT. 



(14) 



The integral can be performed in the complex fre- 
quency plane with the substitution cj — *■ -I- iuj2. It 
is then written 



d-{T)e 



lUJlT -LO2T 



dr. 



(15) 



The factor e^'^^'^ in this integral is bounded in the 
upper half of the complex plane for r > and in the 
lower half plane for r < 0. Since a{T < 0) = 0, this 
means that a{Lu) is analytic for the upper half of the 
complex plane. This means that one can use Cauchy's 
theorem and therefore the relation 



duj' = 



(16) 



holds in the usual way in the upper half of the complex 
plane. (Please see Wootcn (^J, Jackson (I14p or Arfkcn 
(|13P for further details on the integration and derivation 
of the Kramers-Kronig relations). One gets the expres- 
sion 



P 



(17) 



where P as usual denotes the principal part. 

From this the Kramers-Kronig relations can be in- 
ferred. Using the fact that in the time-domain Im 



[^(t)] = 0, which means that a{ 
write a simpler form of them as 



ai =. -P 

2UJ 

02 = P 



dijj 



,uj'a2{uj') 



duj' 



, cri(cj') 



,2' 



(j*{lS) we can 

(18) 
(19) 



The Kramers-Kronig relations are tremendously use- 
ful in the analysis and determination of optical spectra. 
When one knows one component of a response function 
for all frequencies, the other component automatically 
follows. As will be discussed below, they are used exten- 
sively in Fourier Transform Infrared Reflectivity (FTIR) 
measurements to determine the complex reflectivity (am- 
plitude and phase) when only measuring reflected power. 

The above discussion was based on some generalized 
frequency dependent conductivity o(uj) and its model in- 
dependent properties. In real materials, one has a zoo 
of different possible contributions to the electromagnetic 
response. For instance in simple metals with periodic 
translational symmetry (i.e. a crystal), one can iden- 
tify a number of absorptions that satisfy the g = con- 
straint. For the schematic band structure shown in Fig. 
[2] (left), one expects a broad feature at finite energy (in 
red) in the optical conductivity (Fig. |2] (right)), which 
comes from the sum over all possible direct interband 
band absorptions (red arrow) in which electrons are pro- 
moted from below Ej? across an energy gap to a higher 
lying band, with zero net momentum change. These are 
so-called vertical transitions. Near Ej? absorptions with 
low but finite lo (green arrow) are only possible if strict 
translational symmetry has been broken by, for instance, 
disorder. Electrons moving in Bloch waves with mean 
free path £, can violate strict momentum conservation 
in optical absorption at momentum scales on the order 
of 27r/^. One can think of this heuristically as smear- 
ing out the band structure on this scale to allow vertical 
transitions within the same band. This gives a peak cen- 
tered at zero frequency (in green) . In the limit of perfect 
translational symmetry this peak would be a delta func- 
tion centered at w = 0. One can also have excitations of 
harmonic waves of the lattice (phonons) if such phonons 
posses a net dipole moment in the unit cell. They appear 
as distinct and frequently very sharp peaks in the optical 
conductivity (in blue). Many different other excitations 
exist as well that can potentially be created via optical 
absorption. 

Having introduced the general idea of optical response 
functions, I now discuss the derivation of them, both clas- 
sically and quantum mechanically. 



B. Classical Treatments: Drude-Lorentz 

Almost the simplest model of charge conduction we 
can conceive of is of a single charge e of mass m, driven 
by an electric field and subject to a viscous damping 
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Frequency 



FIG. 2 (Color) Various difTerent optical absorptions that sat- 
isfy the g = constraint of optical spectroscopy can appear in 
the optical conductivity of simple metals. Near Ef intraband 
absorptions [green]. Interband absorptions [red]. Lattice vi- 
brations (phonons) [blue] . 



force that relaxes momentum on a time scale r. Consider 
the force equation describing this situation 



— ~eE — mx' /t. 



(20) 



If we assume harmonic motion then x ~ xoe"*'^* and 
E = i?oe~"^*. Substituting in for x and E and solving 
for we get 



erSo 1 

m 1 — ILOT 



(21) 



If we then consider an ensemble of such charges with 
density N, and realize that the maximum current density 
is Jo = Ncx'q, we get the relation 



Jo 



Ne^rE, 







1 



1 — iujT 



(22) 



Using the previously defined relation J — aE we find 
the relation for the frequency dependent 'Drude' conduc- 
tivity is 



Ne'^T 1 
m I — it 



Ne^T 1 -t- iujT 
m 1 + Lo'^T^ 



(23) 



Interestingly, this classical model, first conceived of by 
Paul Drude, is actually of great use even in the analysis 
of particles obeying quantum mechanical statistics. As 
discussed below one finds the same functional form to 
leading order in that case as well. 

The Drude model (plotted in Fig. [3| demonstrates a 
number of important features that are found generally 
in response functions. Its limit at ^ is Ne^r/ra is 
well behaved mathematically and equivalent to the DC 
value. As shown in Fig. |3] it has real and imaginary com- 
ponents that differ from each other. At low frequency 
the current's response is in phase with the driving field 
and purely dissipative (real). At intermediate frequen- 
cies (when the driving frequency is of order the scatter- 
ing rate l/r), the real conductivity falls to half its DC 



value and imaginary conductivity peaks and equals the 
real value. At higher frequencies, two things happen. 
Not only does the current begin to lag the driving field, 
but its overall magnitude decreases as the response has 
a harder time keeping up with the driving electric field. 
The conductivity also exhibits distinct power laws in its 
various frequency limits: linear for the imaginary part at 
low uj, 1/w at high w, and for the real part at high 
frequency. 



■ Real Drude conductivity 
Imaginary Drude conductivity 




Frequency 



FIG. 3 (Color) Frequency dependent Drude conductivity 
with scattering rate l/r = 15 on linear scale (top) log-log 
scale (bottom). 



In the limit of zero dissipation 1 /t 
has the limiting form 

I m 



Equation 23 



CT2 



(24) 
(25) 



(Ti is zero everywhere, but at w = 0. As cti is pro- 
portional to dissipation, this demonstrates that a gas 
of coUisionless electrons cannot absorb photons at finite 
frequency. This can be shown directly by demonstrat- 
ing that the Hamiltonian that describes this interaction. 
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commutes with the momentum operator p and hence it 
has no time dependence. With the exception of umklapp 
scattering, electron-electron interactions do not change 
this situation, as such interactions cannot degrade the 
total system momentum. 

The dependence of the Drude conductivity on r, N , 
and m, shows that the optical mass of charge carriers in 
metals can be obtained if for instance the charge density 
is known from other parameters like the Hall coefficient. 
For instance, in Fig. |4]the mass determined by optical 
conductivity is plotted against the linear coefficient in the 
specific heat (which is proportional to the electronic den- 
sity of states and hence the mass). From simple metals 
to exotic heavy fermion materials, it shows a dramatic 
linear dependence over 3 orders of magnitude. 
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FIG. 4 (Color) Specific- heat values 7 vs effective mass m*, 
evaluated from the optical data using spectral-weight argu- 
ments (I23I1 



It is convenient to express the prefactors of the Drude 
conductivity in terms of the so-called plasma frequency 

The plasma frequency is equivalent to 



the frequency of the free longitudinal oscillations of the 
electron gas The conductivity can be expressed 



a(w) 



47r 1 



(26) 



Therefore within the Drude model, the optical re- 
sponse is fully determined by two frequencies: the plasma 
frequency LOp and the scattering rate l/r. Since ujp is 
many orders of magnitude greater than l/r, this allows 
us to define three distinctly different regimes for the op- 
tical response. 

At low frequencies w <C 1/t, in the so-called Hagen — 
Rubens regime, the conductivity is almost purely real, 
frequency independent, and approximately equal to the 
DC value. At intermediate frequencies, in the relaxation 
regime where u is on the order of 1 /r, one must explicitly 
take into account the lot factor in the denominator of 
Eq. 23 As mentioned above when w = l/r the real 
conductivity falls to half its DC value and the imaginary 
conductivity peaks and is equal to the real value ^. The 
significance of the high frequency regime uj > ujp can 
be seen in rewriting the conductivity as the dielectric 
function. It is 



e(^) = 1 



uj^ — iuj/r 



(27) 



The real and imaginary parts of this expression are 
given in Fig. [5j One sees that for the conventional 
case where tOp ^ l/r, the plasma frequency is the fre- 
quency at which the real part of the dielectric function 
changes sign from negative to positive i.e. Up sets the 
scale for the zero crossing of ei. An analysis of the reflec- 
tion and transmission using the Fresnel equations (see 
Jackson (fT¥)l ), shows that above the zero crossing, met- 
als described by the Drude model become transparent. 
Hence this high frequency limit is called the transparent 
regime. 

Despite its classical nature, the Drude model describes 
the gross features of many metals quite well at low fre- 
quencies. We give a number of examples of its use and 
extensions below. Of course, our interest in the electro- 
dynamics of solids extends far beyond the case of simple 
metals. We will find the Drude model lacking in this 
regard. And even for 'simple' metals, one has the inter- 
esting aspects of finite frequency absorptions that param- 
eterize band structure and, for instance, give copper and 
gold their beautiful colors. The Drude model is wholly 
inadequate to describe such finite frequency absorptions. 
For such processes in more complicated metals, semicon- 
ductors, and insulators, one can model them quantum 
mechanically in a number of ways, but we can also gain 
predictive power and intuition from an extension to the 
classical Drude model called Drude-Lorentz. 

Here we envision the electrons are additionally also 
subject to a simple harmonic restoring force —Kx. Of 
course, such a model is physically inadequate to describe 



* In actual materials, high energy excitations give a contribution to 
the dielectric function that renormalizes the observed plasma 

frequency downward to the screened plasma frequency 

This is the actual frequency of plasmon oscillations in a real 

material. 



Of course in real materials, there are interband excitations which 
can also contribute to the conductivity. The Drude model as 
presented here is an idealized case. 
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FIG. 5 (Color) Frequency dependent Drude dielectric func- 
tion with scattering rate 1/r — 15 and plasma frequency 
Ljp — 60. The real part of the dielectric function changes 
sign at ujp. 





\ Real Drude-Lorentz 

\ conductivity 

\ Imaginary Drude-Lorentz 

\ conductivity 

\ T = 1/15 
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FIG. 6 (Color) Frequency dependent Drude-Lorentz complex 
conductivity with scattering rate 1/r = 15 and oscillator res- 
onant frequency uiq = 45. The real conductivity is zero at 
LO = and positive for all frequencies. 



finite frequency absorption in most band semiconductors 
and insulators, as their insulating nature is due to the 
properties of filled bands and not electron localization. 
Nevertheless, numerous aspects of absorption at 'band- 
edges' in semiconductors can be modeled phenomenolog- 
ically with Drude-Lorentz. Of course its applicability to 
the modeling of harmonic phonon absorptions is obvious. 
With u!q = K/m, we extend Eq. 



20 



as 



mx" = —eE — mx' /t 



(28) 



Proceeding in exactly the same fashion as for the sim- 
ple Drude model, we obtain for the Drude-Lorentz con- 
ductivity 



(29) 



C. The Quantum Case 

The above Drude and Drude-Lorentz models are clas- 
sical models used to describe intrinsically quantum me- 
chanical phenomena and the fact that they are useful at 
all to describe solids is surprising. It is clear that a quan- 
tum mechanical treatment is desired. The most com- 
monly used method for the quantum mechanical calcula- 
tion of electromagnetic response of materials is the Kubo 
formalism. It is based on the fiuctuation-dissipation the- 
orem, which relates the spontaneous fluctuations of a sys- 
tem described via its correlation functions to its driven 
linear response. The derivation here loosely follows the 
one of Wooten ( 1 ) and Dressel and Griiner . 

We calculate conductivity, by calculating the absorp- 
tion rate of a system based on the power dissipated 
P — aiE^ . We first use Fermi's Golden rule to calculate 
the scattering probablility that incident radiation excites 
an electron from one state \s) to another \s') 



One can see that in the limit ujq ^ Q the Drude relation 
is obtained. 

As mentioned above, despite its classical nature and 
the inapplicability of the underlying physical picture for 
some situations, in some cases the Drude-Lorentz model 
can be used to quantify finite frequency absorptions. In 
the case of the band edge absorption in semiconduc- 
tors, a number of oscillators of the form of Eq. |29] 
with different weights can be used to capture the de- 
cidely non-Lorentzian line shape. One can convince one- 
self that since any lineshape can be fit arbitrarily well 
using an arbitrarily large number of oscillators and more- 
over since the imaginary part of the response follows from 
a Kramers-Kroning transform of the real part, it is al- 
ways perfectly feasible to parametrize the response using 
Drude-Lorentz . 



2tt 



\{s'\H,\s)\H{uj-uJs'+u:s). 



(30) 



\s) and \s') are the system's eigenstates, which may 
include the effects of many body physics. The interaction 
Hamiltonian for charge with the electromagnetic field to 
first-order in field is 



e 

2mc 



N 



N 



J2[P^ ■ Ain) + Ain) ■ p,] -eY^Hn), (31) 



i=l 



where p, A, and <I> are the quantum mechanical opera- 
tors for momentum, vector potential, and scalar poten- 
tial and c is the speed of light. Since the vector potential 
A depends on position, it does not in general commute 
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with p. However since p = —ih\/, one can show that 
p ■ A — A ■ p = —ihS7 ■ A and since within the Coulomb 
gauge V ■ A — 0, p and A can in fact commute. This 
means that for purely transverse waves with $ = the 
interaction Hamiltonian simplifies to 



e 

mc 



N 



(32) 



By substituting in for the canonical momentum p — 
mv~eA/c, dropping all terms that arc higher than linear 
in A, and replacing the summations by an integral one 
Hi takes the form 



H^ = -^-JdrJ^{r)-A^{r). (33) 
This gives a Fourier transformed quantity of 



H 



-J^{q)-A^{q). 



(34) 



This means we can substitute in for the form of in 
the expression for Ws~^s' to get 



Ws. 



-I 



s'\Hl\s)Y5{uj-Us' +ujs)- (35) 



The matrix element in this expression follows from the 
form of HT ■ It is 



{s'\Hl\s) 
which leads to 



27r 



{sV^{q)\s){s\J^*{q)\s'), 
\A^{q)\^5{uj -UJs- +LOs)- 



(37) 



where J^*{q) = j'^{—q). One then sums over all occu- 
pied initial and all empty final states W = J^s s' ^s~*s'- 
The dissipated power per unit time and volume at a par- 
ticular photon frequency u follows after a few mathemat- 
ical steps as 



\A^(.q)\'Y: 



UJ 



p = huW 



dt{s\j'^{q,0)J^*{q,t)\s)e~ 



(38) 



The electric field is related to the vector potential as 
E'^ — iujA^ /c for transverse waves. Upon substitution 
we get the absorbed power per unit volume per unit time 
expressed as a current-current correlation function 



P-|£;^(g)P^^ j dt{s\J^{q,Q)J^*{q,t)\s)e-^-\2>9) 

Using our previously given relation P — aiE"^, the 
expression for the real part of the conductivity is 



dt{s\J^{q,Q)J^*{q,t)\ 



s)e 



.(40) 



The imaginary part of cr(a;) follows from the Kramers- 
Kronig relation. This is the Kubo formula. It has the 
form expected for the fluctuation-dissipation theorem, 
that of a correlation function (here current-current) av- 
eraged over all the states \s) and describes fluctuations 
of the current in the ground state. The conductivity de- 
pends on the time correlation between current operators 
integrated over all times. 

The above formalism is general and applies to any set 
of system states \s). An extension can be made in cases 
where Fermi statistics applies. The end result is (see 
Dressel and Griiner ^ for the full derivation) 



m?uj (27r)'^ 



\{s'\p\s)\'Ds>s{nw). 



(41) 



Here (s'|p|s) is known as the dipole matrix element. 
Dgis is the so-called joint density of states, defined as 
^■^Y J S{huj~hujs's)dk. This equation, which is often re- 
ferred to as the Kubo-Greenwood formula, is extremely 
useful for the calculation of higher lying interband tran- 
sitions in metals and semiconductors when the states \s) 
(36) and \s') belong to different bands (as for instance illus- 
trated by the red arrow in Fig. |2|. It also has the form 
that one may naively expect; it is proportional to the 
joint density of states and a matrix element which incor- 
porates aspects like symmetry constraints and selection 
rules for allowed transitions. 



We can use the Kubo formalism (Eq. 40 1 to calculate 
the leading order metallic conductivity in a slightly more 
rigorous fashion than was done in the classical approach 
above. We start with Eq. |40]and then posit as an ansatz 
that finite currents relax as 



J{q,t) = J{q,0)e-'/\ 



(42) 



and assume that the current correlation time r has no 
dependence on q. Inserting this into the Kubo formula, 
one gets that 

^i(9,^) ^^T.J rfte— *-l*l/-(s|j2(g)|5). (43) 

One now needs an expression for the current fluctu- 
ations at flnite q. The expectation value for J is zero. 
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but fluctuations lead to a flnite value of (J^). For 
small q, one can show in the 'dipole approximation' that 
J{q) = — ^ Si Pi where i labels the individual charges. 
With the reinsertion of a complete set of states s' one 
gets the expression 



<Ti(g,w) 



dte-^^t-\t\/r J2 (44) 

S.S' 



The quantity '2J2s,s',^\ ^^^t['l! in Eq. 

called the oscillator strength fs,s'- Here uj in Eq~44| 
is the energy difference between states \s) and |s') as is 
ujs^s' in the definition of the oscillator strength. Below, 
we pay special attention to the oscillator strength when 
we calculate the optical sum rules. 

The oscillator strength can be easily calculated for the 
case of free electrons where the energy is hw = and 

2 2 

the average momentum is |(s'|pi|s)p = ^^j-- In such a 
case /s,s' — N the total number of free charge carriers. 
After performing the integral in Eq. |44]the final result for 
the leading order conductivity calculated via the Kubo 
formalism is 



^Kubo 



1 



(45) 



which is the exact same result as calculated classically 
via the Drude model above. For more details see the 
chapters in Wooten (? ) and Mahan ([3]). 

III. TECHNIQUES 

As mentioned above, the incredibly large spectral 
range spanned by the typical energy scales of solids means 
that many different techniques must be used. Below I 
discuss the most used techniques of microwave measure- 
ments, THz spectroscopies, optical reflectivity, and cUip- 
somctry that are used to span almost six orders of mag- 
nitude in measurement frequency. These various regimes 
are loosely defined below, but of course there are overlaps 
between them ^. 



A. Microwaves 

Until recently measurements in the microwave range 
(100 MHz - 100 GHz) were performed almost entirely 



using microwave cavity resonator techniques. ^ In mi- 
crowave cavity resonance techniques, one is less sensitive 
to such considerations as the sample forms part of a res- 
onance circuit, which dominates the measurement con- 
figuration. Typically one measures transmission through 
a cavity, which is only possible when at resonance. The 
technique is widely used for the study of dielectric and 
magnetic properties of materials in the GHz range, be- 
cause of its high sensitivity and relative simplicity. 

A microwave resonator is typically an enclosed hol- 
low space of cylindrical or rectangular shape machined 
from high-conductivity metal (typically copper or su- 
perconducting metals) with interior dimensions compa- 
rable to the free-space wavelength. On resonance, an 
electromagnetic standing wave pattern is set up inside 
the cavity. One characterizes the resonance characteris- 
tics of the cavity, both with and without a sample in- 
side. Upon sample introduction the resonance's center 
frequency shifts and it broadens. From the shifting and 
broadening, one can define a complex frequency shift, 
which with knowledge of the sample size and shape al- 
lows one to quantify the complex conductivity directly. 
For instance, in the limit of very thin films, the center 
frequency shift is proportional to cr2 and the broadening 
is proportional to ai . For samples much thicker than the 
skin or penetration depth, different relations apply ()25p . 
The fact that the measurements are performed "on res- 
onance" means that a very high sensitivity is achieved. 
The technique does have the considerable disadvantage 
that only discrete frequencies can be accesses, as one is 
limited to the standing wave resonance frequencies of the 
cavity, which means that true spectroscopy is not possi- 
ble. In such measurements one typically determines a 
sample's complex conductivity as a function of tempera- 
ture at some flnite flxed frequency (See Figs. 32 33 and 
34|. 



In contrast , the Corbino geometry is a measurement 
conflguration that is capable of broadband microwave 
spectroscopy. It is compatible with low T and high fleld 
cryogenic environments and capable of broadband mi- 
crowave spectroscopy measurements from 10 MHz to 40 
GHz. It is a relatively new technique, but one which 
has recently been used with great success by a number 
of groups in the correlated electron physics community 
(Uni [571 |5H1 [ini [Sni E]) • Although powerful, it has the dis- 
advantage that its use is conflned to samples which have 
typical 2D resistance within a few orders of magnitude 
of 50 Ohms. 



Although I don't discuss them, one might naturally include vari- 
ous radio frequency techniques in this list (which would increase 
the frequency range by another three orders of magnitude at 
least at the low end). At typical experimental temperatures, 
the low frequency available in techniques like mutual inductance 
measurements allows almost the DC limit of quantities like the 
superfluid density in superconductors to be measured. Interested 
readers are directed to Ref. (I24II and reference therein. 



^ At lower frequencies (kHz and MHz), measurements can be per- 
formed by attaching contacts to samples and the complex con- 
ductivity can be measured by lock-in techniques, network analyz- 
ers, and impedance analyzers. Such methods become problem- 
atic in the GHz range however, as wavelengths become compara- 
ble to various measurement dimensions (sample size, microwave 
connector dimensions, cable lengths), and capacitive and induc- 
tive effects become appreciable. 
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FIG. 7 (Color) View of the UCLA He3 'CryoMag' system. 
Cavities can be interchanged quickly to enable a number of 
different frequencies to be accessed. 



In a Corbino geometry spectrometer a microwave sig- 
nal from a network analyzer is fed into a coaxial trans- 
mission line. A schematic is shown in Fig. |8] (left). The 
signal propagates down the coaxial line and is reflected 
from a sample that terminates an open ended coaxial 
connector. The network analyzer determines the com- 
plex reflection coefficient of the sample which can be 
related to the complex sample impedance Zi^ via the rela- 
tion Sii = where Zq is the coaxial line impedance 
(nominally 50 fi). Because the sample geometry is well 
defined, knowledge of the sample impedance yields in- 
trinsic quantities like the complex conductivity. In the 
typical case of a thin film where the skin or penetration 
depth is much larger than the sample thickness d, the 
sample impedance is related in a straightforward fashion 

to the complex conductivity as cr = ^'^^'^^^^^'^ where ri 
and r2 are the inner and outer conductor radii respec- 
tively. 

An experimental challenge in Corbino geometry mea- 
surements is that the coaxial cables and other parts of 
the transmission lines can have strongly temperature de- 
pendent transmission characteristics. Errors in the in- 
trinsic refiection coefficient, coming from standing wave 
refiections or phase shifts and damping in the transmis- 
sion lines, are accounted for by performing a number of 
calibration measurements. It is imperative that the same 
cryogenics conditions are reproduced between calibration 
and each subsequent measurement. In general this is eas- 
ily done by ensuring the same starting conditions and 
using a computer controlled cool-down cycle. Recently 
it has been demonstrated that a three sample calibra- 
tion (open, short, and a standard 50 ohm resistor) at all 
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FIG. 8 (Color) (left) Experimental geometry for a Corbino 
spectrometer, (right) The real and imaginary parts ai+ia2 of 
the optical conductivity spectrum of UPd2 AI3 at temperature 
T — 2.75 K. The data are fit to <Jdc{^ — iior)'^ with Odc = 
1.05X 10^(n — cm) ^ and r = 4.8 x 10 ^^s and show excellent 
agreement with the Drude prediction. Adapted from Refs 
(IMllZl) 



temperatures dramatically increases the precision of the 
technique over previous single sample calibrations ((26|l . 

As an example of the power of the technique we show 
in Fig. [s] (right) recent measurements taken on films of 
the heavy-fermion compound UPd2Al3 at T = 2.75 K. 
The data show almost ideal Drude behavior with a real 
conductivity of Lorentzian shape and an imaginary con- 
ductivity that peaks at the same uj as the real conductiv- 
ity has decreased by half. The anomalous aspect is the 
Drude peak is seen to be of remarkably narrow width - 
almost 500 times narrower than that of a good metal like 
copper. This is a consequence of the mass renormaliza- 
tions in the heavy fermion compounds and a feature that 
would have been completely undetectable with a conven- 
tional spectrometer. 

Recently there also been the very interesting develop- 
ment of a bolometric technique for high-resolution broad- 
band microwave measurements of ultra-low-loss samples, 
like superconductors and good metals (j32|l . This tech- 
nique is a non-resonant one where the sample itself is 
actually used as a detector in a bolometric fashion. Mi- 
crowaves are fed through a rectangular coaxial transmis- 
sion line into which the sample under test and a refer- 
ence sample are mounted. Small changes in the sample 
and reference's temperature are monitored as they ab- 
sorb microwave radiation. The key to the success of this 
technique is the in situ use of this normal metal ref- 
erence sample which calibrates the absolute microwave 
incident power. As the absorbed power is proportional 
to the surface resistance Rs {lo) , the independent measure 
of the incident power allows one to measure the surface 
resistance precisely. 

In such an apparatus, the sample temperature can be 
controlled independently of the 1.2 K liquid-helium bath, 
allowing for measurements of the temperature evolution 
of the absorption. The minimum detectable power of 
this method at 1.3 K is 1.5 pW, which corresponds to 
a surface impedance sensitivity of ~ Ififl for a typical 
1mm X 1mm platelet sample. The technique allows very 
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sensitive measurements of the microwave surface resis- 
tance over a continuous frequency range on highly con- 
ducting samples. This is a region of sample impedances 
generally inaccesible with the Corbino technique and only 
available at discrete frequencies in microwave cavities. 

A disadvantage of the technique is that although one 
measures i?s(^) the quantity of interest is typically the 
complex conductivity <7{i-u). One must make various 
assumptions to get this from the relation for surface 
impedance 



Zs = Rs + iXs = 



(Tl + i(T2 



(46) 



as the surface reactance Xs is not measured. At low tem- 
peratures and frequencies in the superconducting state CT2 
is mostly determined by the superconducting response 
and can be related to the independently measured pen- 
etration depth A. a I then follows from the relation 
Rs = ^fiQiu^X^cTi. At higher temperatures and frequen- 
cies a more complicated iterative procedure must be used 
that is based on the Kramers-Kronig transform (f5^ . 

This technique has been used to measure the real part 
of the conductivity of YBa2Cu306.5 single crystals across 
a wide frequency range as shown in Fig. [9] One can see a 
quasi-Drude like peak at low frequencies that comes from 
quasi-particle excitation at the nodes of the d-wave su- 
perconductor. The peak has a narrow width of the order 
of 5 GHz, signifying a collapse of quasiparticle scattering 
in the superconducting state and very long quasiparti- 
cle scattering lengths. This is indicative of the very high 
quality of these YBa2Cu306.5 single crystals. Measure- 
ments of a broadband nature on these kind of materials 
are not possible with other techniques. As powerful as 
it is, the technique is hampered by the fact that various 
assumption must be made about a2- This is problematic 
for materials where its dependence is not known a priori. 

In addition to the ones detailed above, there are 
many other techniques which have been used for ma- 
terials characterization in the microwave range. These 
include microstrip and striplinc resonators, transmission 
line wave bridges, radio frequency bridges, parallel-plate 
resonators, and hehcal resonators. See Ref. ([21 |33l I112p 
for further details. 



B. THz spectroscopy 

Time scales in the picosecond (10""'^^ sec) range are 
among the most ubiquitous in condensed matter systems. 
For example, the resonant period of electrons in semicon- 
ductors and their nanostructures, the scattering times of 
electrons in metals, vibrational frequencies of molecular 
crystals, superconducting gap energies, the lifetime of bi- 
ologically important collective vibrations of proteins, and 
- now - even the transit time for an electron in Intel's 
new THz transistor - these are all picosecond phenom- 
ena. This ubiquity means that experimental probes em- 
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FIG. 9 The real part of the microwave conductivity cri{u}) 
extracted from measurements of Rs{(jj) as described in the 
text and Ref. ((32)l . Figure adapted from l(32)l . 



ploying Terahertz (THz) electromagnetic radiation are 
potentially quite powerful. It is unfortunate then that 
this spectral range lies in the so-called 'Terahertz Gap' 
- above the capabilities of traditional electronics and the 
microwave range, but below that of typical optical in- 
strumentation. The THz spectral region has been a tra- 
ditionally difficult part of the electromagnetic spectrum 
to work in because of a number of reasons including weak 
sources, long wavelengths, and contamination by ambient 
room temperature black body radiation. In recent years, 
however, there have been a number of developments that 
allow measurements in the THz range in a manner that 
was not previously possible. These come in the form of 
time domain THz spectroscopy and the increasing preva- 
lence of Backward Wave Oscillator (BWO) based spec- 
trometers. In recent years, however, a number of dra- 
matic technical advances such as time-domain THz spec- 
troscopy (TDTS) using so-called 'Auston' switch genera- 
tors and detectors have enabled measurements that span 
this gap in measurement possibilitites. 

THz spectroscopy has become a tremendous growth 
field (|34[l . finding use in a multitude of areas including 
characterization for novel solid-state materials ((551 
optimization of the electromagnetic response of new coat- 
ings (P7|) . probes of superconductor properties (j551 
security applications for explosives and biohazard detec- 
tion (|40)l . detection of protein conformational changes 
(I4ip . and non-invasive structural and medical imaging 



TDTS works (See Ref. (ji^ for an additional excellent 
short summary) by the excitation of a source and acti- 
vation of a detector by ultrafast femtosecond laser (typ- 
ically Ti-sapphire) pulses. A basic schematic is shown 
in Fig. |10[ A femtosecond laser has, via a beam split- 
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FIG. 10 A diagram of a typical experimental layout used in 
TDTS. 
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FIG. 12 (Color) (left) Time-domain trace of electric field 
transmission through 60 pm thick SrTiOa subtrate. Multiple 
reflections from front and back surfaces are readily visible, 
(right) Power spectrum of the transmission through the same 
60 pm thick SrTiOa subtrate. The data is obtained from the 
squared Fourier transform of the data in (left) and ratioed to 
a reference aperature signal. From Ref. (|46j) 




FIG. 11 An optical image of the two-contact photoconductive 
'Auston' switch antenna structure used in TDTS. Two elec- 
trodes are grown on top of an insulating high defect semicon- 
ductor, like low temperature grown GaAs or radiation dam- 
aged silicon. Electrode spacing is typically on the order of 10 
microns. 



ter, its radiation split off to fall on source and detec- 
tor Auston switches, which are typically pieces of low 
temperature-GaAs or radiation damaged silicon with two 
electrodes grown on top in a dipole arrangement and 
separated by approximately 20 /im (Fig. 111. Before 
laser illumination, the switch has a resistance of a few 



megaohms. After the fast illumination by the femtosec- 
ond laser pulse, the source switch's resistance falls to a 
few hundred ohms and with a bias by a few tens of volts, 
charge carriers are accelerated across the gap on a time 
scale of a few picoseconds. Their acceleration produces a 
pulse of almost single-cycle radiation, which then propa- 
gates through free space and - coUimated by mirrors and 
lenses - interacts with the sample. Measurements have 
been typically performed in transmission. After passing 
through the sample, THz radiation falls on the second 
Auston switch, which is also activated by the femtosec- 
ond pulse. Whereas the first Auston switch was DC bi- 
ased, the 2nd Auston switch is biased by the electric field 
of the THz radiation falling upon it. Current can flow 
only across the second switch with a direction and mag- 
nitude proportional to the transient electric field at the 
instant the short femtosecond pulse impinges on it. The 
experimental signal is proportional to this current with 



a magnitude and polarity that reflects the THz electric 
fleld at the switch. A delay line is then advanced and 
the THz electric field at a different relative times can be 
measured. In this way, the entire electric field profile as 
a function of time can be mapped out as shown in Fig. 
p^lcft). The time-domain pulse is then Fourier trans- 
formed to get the complex electric field as a function of 
frequency. A reference measurement is also performed 
on an aperture with no sample. The measured transmis- 
sion function is the ratio of the signal transmission to the 
reference [Fig. 12 'right)]. 

There are a number of unique aspects to TDTS that 
allow it to work exceptionally. Since the detected signal 
is proportional to the instantaneous electric field, and 
not the power, the measured transmission function is the 
complex transmission coefficient for the electric field, 
which has both an amplitude and a phase. This allows 
one to invert the data directly to get the real and imagi- 
nary optical constants (e.g. the complex conductivity) of 
the material. This is essential for a thorough characteri- 
zation of material's properties. Typical optical measure- 
ments discussed below measure only reflected or trans- 
mitted power, quantities in which the complex optical 
constants of interest are confolded into in a non triv- 
ial fashion. In such measurements one has to measure 
over the largest frequency range possible (even if one is 
only interested in a limited spectral range), extrapolate 
to both DC and infinite frequency, and then Kramers- 
Kronig transform to get phase information. In TDTS 
one can get both components by direct inversion. 

TDTS is also capable of unprecedentedly high signal- 
to- noise ratios. Typically efforts in the THz and far in- 
frared spectral range are complicated by a very large 
ambient black body radiation background. In the case 
of TDTS the fact that the signal is 'time-gated' and co- 
herent, whereas black-body radiation is incoherent, al- 
lows a very high detection efficiency. Additionally, since 
the detected quantity is actually electric field and not 
power (proportional to electric field squared), the noise 
in the spectral power is reduced by a factor proportional 
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FIG. 13 The dynamic (frequency dependent) phase stiffness 
Te{uj) as a function of temperature T. Data are shown for 
two samples, one with Tc 33 K. The data for each sample 
correspond to measurement frequencies of 100, 200 and 600 
GHz. Each set of curves identify a crossover from frequency- 
independent to frequency-dependent phase stiffness. The 
dashed line shows that the crossover is set by the Kosterlitz- 
Thouless-Berenzkii condition for two-dimensional melting. 
From Ref. (|39} 



to the electric field itself. These aspects allow detection 
of transmission signals approaching one part in 10^ . This 
is an incredibly large dynamic range and is essentially un- 
precedented in optically based spectroscopies where one 
part in 10'^ is typically considered extremely good. 

TDTS was used for instance in the measurement of 
thin Bi2212 cuprate superconducting films to show ev- 
idence for phase fiuctuations above Tc. Corson et al. 
measured transmission through thin films and using the 
unique phase sensitivity of TDTS could show that there 
was an enhanced (T2 (Fig- 13 1 even at temperatures well 
above Tc . This was interpreted as a finite phase stiff- 
ness on short length and time scales. Superconducting 
fluctuations as such should manifest themselves as a fl- 
nite superfluid stiffness measured at flnite frequency. The 
imaginary part of the optical conductivity is sensitive to 
superfluid stiffness through the relation a2 — cq ^j^J" , 

where ctq = ^ is the quantum of conductance for 
Cooper pairs divided by a length scale (typically the c- 
axis lattice constant) and Tg is the generalized frequency 
dependent superfluid stiffness. Here, the superfluid stiff- 
ness Tg is an energy scale expressed in temperature units. 
The stiffness is the energy scale to introduce phase twists 
in the superfluid order parameter \I' = Ae"^. 

The other major development for THz spectroscopy 
is the increased prominence of Backward Wave Oscilla- 
tors (BWOs) Although developed in the 60's, BWOs are 
gaining increasing prominence in the investigation of cor- 
related systems in the important THz range (|47| . These 
BWOs are traveling wave guide tubes, capable of pro- 



ducing very monochromatic THz range radiation over a 
relatively broad range per device. With a number of 
tubes, it is powerful method to cover the 0.03 THz to 1.5 
THz range in a continuous wave conflguration. 

The longer wavelengths in this spectral region give the 
capabilities to measure phase of transmitted waveforms 
through Mach-Zender two-beam polarization interferom- 
etry. This method also allows direct precision measure- 
ments of both real and imaginary components of the com- 
plex optical response without resort to Kramers-Kronig 
transforms and their associated ambiguities. A spectrom- 
eter based on BWOs is easily integrated with a low tem- 
perature cryostat and magnetic field system as well as 
glove boxes for environmental control. 



C. Infrared 

Fourier transform infrared reflectivity (FTIR) is the 
workhorse spectroscopy for optical characterization of 
solids (Fig. 151. Measurements are possible from ~ 10 
cm~^ to approximately 25,000 cm~^, although measure- 



ment get increasingly difficult below 50 cm~^. Heroic 
efforts can push the lower end of this range slightly be- 
low 10 cm~^ ()48)) . Here the quantity of interest is typ- 
ically the power reflectivity (transmission measurements 
are possible as well). One shines broadband light, some- 
times from several different sources, on a sample surface 
and via an interference technique measures the reflected 
intensity. As the source spectrum and detector response 
can have all kinds of frequency structure, it is necessary 
to reference the reflected signal against a standard sample 
to obtain the absolute reflectivity. Typically the reference 
spectrum is chosen to be a noble metal like gold whose 
reflectance can to good approximation be taken to be 
unity over a broad frequency region below the material's 
plasma frequency (see below). Typically gold is evapo- 
rated onto a sample as shown in Fig. [l4]or the sample is 
replaced with a mirror to perform this referencing. 

As the materials we are interested in typically have 
interesting experimental signatures in both the real and 
imaginary response, it is not sufficient to to compare to 
theory by simply measuring quantities like absorbence. 
We want quantities like the complex conductivity or di- 
electric function. It is not necessarily straightfoward to 
obtain a complex quantity from the measured scalar mag- 
nitude of the reflectivity R{u!). Unlike TDTS, phase is 
not measured in FTIR measurements. Therefore typi- 
cally use is made of the above discussed Kramers-Kronig 
transforms, which apply to any causal response function. 
If one knows one component of the response (real or 
imaginary) for all frequencies then one can determine the 
other component. For reflected power, one can obtain the 
reflected phase as 



(47) 
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Reflectivity (Far-IR - Near-IR) 




FIG. 14 (Color) Absolute reflectivity R{lo) in FTIR is ob- 
tained by measuring the reflected intensity and ratioing with 
the intensity after gold evaporation or from a reference mirror. 
Figure courtesy of A. Kuzmenko. 

Of course the problem is that one doesn't measure over 
an infinite frequency range. In typical FTIR spectoscopy, 
the usual mode of operation is therefore to measure over 
as large a frequency range as possible and then extrap- 
olate with various schemes to w — > and w — *■ oo. This 
method works quite well for some materials, although it 
can generate large errors when, for instance, reflectivi- 
ties approach unity as they do in good conductors at low 
frequencies. It can also be rather inefficient and cum- 
bersome as one must measure out to the multi-eV range 
even if one is interested in meV level phenomena. 

As mentioned, FTIR spectoscopy is a standard spec- 
troscopic technique for the characterization of materials 
and chemicals. It is used routinely to identify the pres- 
ence of various chemical bonds in chemistry as in, for 
instance, C = C or C = O, which all have distinct fre- 
quencies and absorptive strengths. It is also the standard 
optical tool for probing materials like high-Tc supercon- 
ductors @ . It has the advantage of being relatively easy 
to perform and possessing a very large spectra range. It 
has the above mentioned disadvantage of only measur- 
ing a scalar - the power - and having a lower bound on 
the spectral range that is at the limit of that explicitly 
relevant for many correlated materials. 



D. Visible and Ultraviolet 

The technique of optical reflection as outlined above 
begins to become more difficult at frequencies on the or- 
der of the plasma frequencies of metals. To do reflectiv- 
ity measurements, one always needs a standard sample. 
Standard metal references cannot be used above their 
plasma frequencies as their reflectivity changes quickly 
with frequency. Quartz or silicon which have an approx- 
imately flat (but low) reflectivity can be used as a ref- 
erence in a grating spectrometer (|49p . However, more 
common nowadays is to use techniques such as ellipsome- 




FIG. 15 (Color) Modified Fourier Transform Spectrometer 
Bruker 66. As shown, the spectrometer is setup to perform 
reflection, but with additional mirrors near-normal incidence 
reflection is possibe. The apparatus works by shining broad 
band light from a variety of sources on a sample. Reflected 
or transmitted intensity is measured with a few different de- 
tectors. A beamsplitter/interferometer apparatus allows the 
resolution of distinct frequency components. Figure from 
http : //www .pi 1 .physik.uni-stuttgart . de/research 
/Methoden/FTIR_e .php. Figure courtesy of N. Drichko. 



try to determine complex optical coefficients in this spec- 
tral range as one can measure complex response functions 
directly. 

Ellipsometry is an optical technique for determining 
properties of surfaces and thin films. Although the 
method was originally used as far back as 1887 by Paul 
Drude to determine the dielectric function of various met- 
als and dielectrics, it only gained regular use as a charac- 
terization tool in the 1970s (50 ; 51). It is now widely used 
in the near infrared (NIR) through ultraviolet (UV) fre- 
quency ranges in semiconductor processing for dielectric 
and thickness characterization. 

When linearly polarized light is reflected from a sur- 
face at glancing incidence the in- (p) and out-of-plane 
(s) of incidence light is reflected at different intensities 
as well as suffering a relative phase shift. The reflected 
light becomes elliptically polarized as shown in Fig. [16] 
The shape and orientation of the ellipse depend on the 
angle of incidence, the initial polarization direction, and 
of course the reflection properties of the surface. An el- 
lipsometer measures the change in the light's polariza- 
tion state and characterizes the complex ratio p of the 
in- (rp) and out-of-plane (tj) of incidence reflectivities. 
The Fresnel equations allow this quantity to be directly 
related to various intrinsic material parameters, such as 
complex dielectric constant or layer thicknesses. 

The ellipsometer itself is designed to measure the 
change in polarization state of the light reflected from 
a surface at glancing incidence. From a knowledge of 
the orientation and polarization direction of the incident 
light one can calculate the relative phase difference. A, 
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two orthogonal directions, as can be seen by construction 
in Fig. 17 In order to get the phase information, the el- 



^sample 



plane^f incidence 



FIG. 16 (Color) Schematic showing the basic principle of el- 
lipsometry. Glancing incidence light linearly polarized with 
both s and p components is reflected with different intensities 
and a relative phase shift. A characterization of the resulting 
elliptically polarized reflected light's minor and major axes, 
as well as tilt angle gives a unique contribution of the complex 
optical constants of the material. 



and the relative amplitude difference (tan^*) between the 
two polarization components that are introduced by re- 
flection from the surface. 

Given the complex ratio p = rp/r^ = e^^tan^ of the 
in- (rp) and out-of-plane reflectivities (rg), various intrin- 
sic material parameters can be generated via the Fresnel 
equations. In a homogeneous sample the complex dielec- 
tric constant e = ei + ie2 is related to p as given in Eqs. 
48 49 and 50 below. Here as usual ei parameterizes the 



polarizability of a material, whereas 62 parameterizes the 
dissipative properties, ez is the complex dielectric con- 
stant perpendicular to the reflection surface and e^. is the 
complex dielectric constant in the plane of the reflection 
surface, and 6 is the angle of incidence as shown in Fig. 

m 



p = Vp/rs = tan^fe*^. 



tz ^sin^6' — ^/e7cos0 



\/\ — ez ^sin^0 + Y^cosfl 
COS0 — A/e,, — %\T?Q 



cos6 



(48) 
(49) 

(50) 



In a typical conflguration, monochromatic light is in- 
cident on a glancing trajectory 6 close to the Brewster 
angle (usually ~ 65 — 85° for bad metals at high fre- 
quencies) with a linear polarization state of 45°. As the 
typical detector measures power and not electric fleld, it 
is essential that the light's elliptical polarization is mea- 
sured over at least 180° (and more typically 360°). To 
completely characterize the phase, it is clearly not suffi- 
cient to simply measure the projection of the ellipse along 



lipse's orientation and major and minor axes must be 
known. Generically, the light originally linearly polar- 
ized at 45° is changed to an elliptical polarization with 
its major axis displaced from 45°. A characterization of 
this ellipse's major and minor axes, as well as tilt angle is 
a direct measure of the complex amplitude reflection co- 
efficients, which can then be related to intrinsic material 
quantities. 




FIG. 17 (Color) To completely characterize the electric fleld 
vectors and their phase, it is insufficient to measure the pro- 
jection of the ellipse along two orthogonal directions. The 
complete ellipse must be mapped out to get its major and 
minor axes and tilt angle. 

The fact that ellipsometry measures the ratio of two 
simultaneously measured values gives it several distinct 
advantages as a characterization tool over simple reflec- 
tivity. It is highly accurate and reproducible even at 
low intensities as many systematic errors are divided 
out. Moreover, unlike reflectivity measurements, the 
technique is self-normalizing and a reference sample is 
not necessary. This is important as the choice of a ref- 
erence can be problematic in the optical and ultraviolet 
spectral range or when the surface is of unknown quality. 
The method is also particularly insusceptible to source 
intensity fluctuations as they are also divided out. Per- 
haps most important for the overall utility of the tech- 
nique is that the method also measures a phase, which 
gives additional information regarding materials proper- 
ties. The phase information can be used to generate the 
complex optical constants (the complex dielectric con- 
stant for instance, or the complex index of refraction n 
and k or the complex conductivity a) which is essential 
information for a complete characterization of a mate- 
rial's optical response. If the optical constants are known, 
the phase information can be used to sensitively measure 
films thicknesses. If conventional ellipsometry has a dis- 
advantage it is that with typical grating monochromator 
data acquisition (i.e. one frequency stepped at a time) it 
is relatively slow as both the grating and polarizers must 
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continue to be moved. In such a setup one loses the rapid 
acquisition time of the multiple frequency multiplexing 
in, for instance, Fourier transform infrared reflectivity 
(FTIR) ®. As mentioned above, the technique is well es- 
tablished and widely used in the NIR through ultraviolet 
(UV) frequency ranges (|52l [53]l . There have been a num- 
ber of attempts to extend ellipsometry to lower infrared 
frequencies (52^ i53i i54,) . Unfortunately the lack of suffi- 
ciently intense sources (among other reasons) has meant 
that such efforts have met with limited success, although 
synchrotron based efforts have had some success in this 
regard ((55|) . 

The energy scales of the visible and ultraviolet, which 
are measured by ellipsometry are not usually explicitly 
relevant for strongly correlated systems. Typically we are 
interested in much lower energies that are on the order of 
the temperatures that phenomena are exhibited. How- 
ever, such energies are relevant for the determination of 
important band structure parameters on, say, 3D mate- 
rials that photoemission cannot be performed on. Exci- 
tations in this range also determine the high frequency 
dielectric constant eoo , which can have a significant (but 
indirect) effect of the low energy physics. Characteri- 
zation of optical constants in the visible and UV range 
is also imperative to constrain the phase information at 
low energies when performing Kramers-Kronig with re- 
flectivity data. Recently an extremely powerful method 
for combining different data sets using a Kramers-Kronig 
consistent variational fitting procedure was developed by 
A. Kuzmenko (j56p . It allows a general procedure for 
combining say, DC resistance, cavity microwave measure- 
ments at a few distinct frequencies, IR power refiectivity, 
and ellipsometry to extract over a broad energy range of 
all the significant frequency and temperature dependent 
optical properties. 



IV. EXAMPLES 
A. Simple Metals 

Noble metals like silver, copper, and gold provide a 
good first step for the understanding the electrodynamics 
of solids. Compare the refiection spectra of silver and 
gold (and aluminum) given in Fig. 18 as a function of 
energy. The refiection spectra of silver shows a sharp 
plasma edge around 3.8 eV Gold has a much lower plasma 
frequency with its plasma edge at a wavelength of around 
2.5 eV. Also shown is the refiection spectra for aluminum 
which has a much larger plasma frequency than any of 
these at around 15.5 eV. The small dip in Al around 1.65 
eV is caused by an interband absorption (see below) and 
is not its plasma frequency. 
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FIG. 18 (Color) Reflectance curves for aluminium (Al), silver 
(Ag), and gold (Au) metal at normal incidence as a function 
of incident energy. The plasma frequencies of silver and gold 
are clearly seen. The small dip in aluminum's curve around 
1.65 eV is caused by an interband absorption and is not its 
plasma frequency. Aluminum's plasma frequency lies at much 
higher energies at around 15.5 eV. Adapted from Ref. (|57p . 



From these spectra we can immediately see the ori- 
gin for the visual difference between silver and gold. As 
mentioned above, the reflection/transmission properties 
of metals greatly change above their plasma frequencies. 
If not for the presence of interband transitions, metals 
would be transparent in this regime. The lower plasma 
frequency of gold is what determines its yellowish color in 
reflection. It reflects less in the blue/ violet portion of the 
spectra and hence looks yellowish. In contrast the plasma 
frequency of silver is in the UV range and so the reflec- 
tion of silver is almost constant at constant throughout 
the infrared to visible range. Aluminum's plasma fre- 
quency is higher yet still and so the major differences in 
its optical properties as compared to silver are invisible 
to the human eye. 

As discussed above, within the Drude model with only 

free carriers the plasma frequency is set by ujp = ' '^'^^'^^ 



This is frequency of the dip in the reflectivity, the zero 
crossing in ei, and the frequency of free longitudinal 
charge oscillations (the plasmon frequency). However, 
the presence of interband excitations, which contribute 
a high frequency dielectric constant eoo changes the sit- 
uation. In real materials, these features are exhibited 
at the so-called screened plasma frequency set by ujp — 

It is interesting to note that the carrier density 



In FIR ellipsometers one typically uses the interferometer of an 
FTIR spectrometer, which retains this multiplexing. 



and masses of silver and gold are almost identical and so 
the difference in their plasma frequencies comes almost 
entirely from the renormalization effects of different in- 
terband transitions eoo. 



Data like that shown in Fig. 18 can, with appropriate 
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FIG. 19 Real and imaginary dielectric functions for silver 
(top) and gold (bottom). The thickness of the curve indicates 
the experimental uncertainty, which originates mainly in the 
Kramers-Kroning transform. From Ref. I|58p . 



extrapolations to high and low frequency, be Kramers- 
Kronig transformed to get the real and imaginary dielec- 
tric function as shown in Fig. 



19 for silver and gold. As 
expected they both show strong almost perfect u) — Q 
Lorentzian Drude peaks of almost equal weight, but gold 
shows stronger interband absorptions, gives it a larger 
that renormalizes its plasma frequency ^ . Such spectra 
can be decomposed into various contributions as shown 



^ Recall from section [V.Bj that a low frequency oscillator (absorp- 
tion) can give a large too ■ AH other things being equal, the lower 
the frequency of the oscillator the higher the eoo. 



in Fig. 20 for silver. As expected there is a strong zero 
frequency Drude peak contribution as well as interband 
contributions. Although the interband piece shows some 
resemblance to the Drude-Lorentz plots in Fig. |6] the 
shape is a more rounded cusp, which is characteristic of 
transitions within the manifold of d-electron states. 




FIG. 20 Decomposition of the dielectric function of silver into 
Drude and interband contributions. The threshold absorption 
for interband excitations is indicated by uji. e^^' is the Drude 



contribution. 5e\ is the interband piece. From Ref. (|60|) . 



In Fig. 



21a similar plot of the complex response func- 
In addition to the metallic 



tions for aluminum is shown 
LD — Q Lorentzian Drude peak, the prominent absorption 
that gives the dip in the reflectivity curve around 1.65 eV 
can be clearly seen. That this dip in the reflectivity does 
not reflect the plasma frequency can be see in the plot 
of the loss function Im {—\), which shows a strong peak 
only at the 'screened' plasma frequency 15.5 eV (i.e. the 
zero cross of ei. 



B. Semiconductors and Band Insulators 

In semiconductors and band insulators, the real part 
of the conductivity is dominated by interband transitions 
(red in Fig. [2]). At low temperatures in a clean undoped 
insulator, naively one expects to see a gap with no con- 
ductivity and then a sudden onset at the gap edge. In 
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FIG. 21 (top) Real and imaginary dielectric functions for 
aluminum, (bottom) Real part of the conductivity and the 
loss function. In addition to the Drude peak, a prominent 
absorption is found around 1.65. This feature gives the dip in 
the reflectivity curve around the same energy. The loss func- 
tion shows a strong peak at the screened plasma frequency 
15.5 eV. From Ref. (f59|l . 



fact, various other contributions are also possible inside 
the gap. Electron-hole pair bound states can form as ex- 
citons and give dissipative response and be seen as sharp 
peaks below the gap. This is found for instance in Fig. 
22 for the ionic insulator KCl, where in addition to the 
band edge excitation there are a series of below gap ex- 
citon absorptions. Even lower measurement frequencies 
than those shown would reveal a series of sharp peaks 
from optical phonons down in the 50 meV range. Such 
behavior can be seen in Fig. [23] for MgO and FeSi. 



If one starts to dope into a typical semiconductor like 
phosphorus into silicon, massive changes in the optical 
spectra occur. At very low doping levels one can discern 
the sharp absorptions of individual donors as seen in the 
bottom of Fig. 
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As the dopant density is increased 
inter-pair and inter-cluster quantum tunneling broadens 
the spectra until a prominent a prominent impurity band 
forms in the gap as shown in Fig. |25j Charge transitions 
can occur in this manifold of impurity states. Below a 
certain critical doping such a doped state retains its in- 
sulating character and is termed an electronic glass. I 
discuss this state of matter below. 
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FIG. 22 Imaginary dielectric function of KCl. The absorption 
around 9 eV is the bad gap. The series of sharp features at 
energies below the band gap are due to excitons. From Ref. 
(ED 



C. Electron glasses 

Many crystalline and amorphous semiconductors show 
a metal-insulator transition (MIT) as a function of 
dopant concentration or stochiometry. On the insulat- 
ing side of the transition, the materials are insulators not 
principally for the reason of completely filled bands (as 
in the case of a band insulator) or interactions (as in the 
case of a Mott insulator), but because electronic states 
at the chemical potential are localized due to disorder. 
Such materials which have a random spatial distribution 
of localized charges have been called - in analogy with 
structural glasses - electron-glasses and are more gener- 
ally termed Anderson insulators. 

At finite temperatures these materials conduct under 
DC bias via thermally activated transport of charges 
hopping between localized states in an impurity band. 
Hence, the DC conductivity tends to zero as the tem- 
perature tends to zero. However, because their insulat- 
ing nature derives from localization of the electronic or- 
bitals and not necessarily a vanishing density of states at 
Ei?, such materials can have appreciable conductivity at 
low frequency deriving from transitions between localized 
states inside the impurity band. Although electron trans- 
port in electronic glasses is still not completely under- 
stood, different regimes for AC transport can be distin- 
guished depending on the strength of electron-electron in- 
teractions, the energy scale being probed and the amount 
of disorder. See Helgren et al. fB5) for a complete dis- 
cussion on this subject. 

Far away from the MIT and in regimes when electron- 
electron interactions are not significant, one considers a 
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FIG. 23 Optical conductivityof MgO (top panel) and FeSi at 
T—iK (bottom panel) showing the strong phonon absorption 
found in such insulating compounds. In the insets are the 
integrated spectral expressed in units of effective charge (See 
Sec. 



V.Ai 



with the electronic contribution removed. For 
MgO the effective measured transverse charge is 1.99, which 
is in good agreement with the formal valences of the Mg"*"^ 
and 0~^ ions. From Ref. I|62[l . For further discussion on 
effective charge see Ref. (|7|). 



so-called Fermi glass. This is an ensemble of localized 
charges whose properties are primarily determined by 
Fermi statistics. When the long-range Coulomb inter- 
action is of the same order as the disorder potentials, 
a so-called Coulomb glass emerges, and finally, near the 
transition to the metallic state, fluctuations lead to a 
quantum critical regime. For these various regimes, cer- 
tain characteristic power law and exponential functional 
forms for the AC and DC transport are expected. For in- 
stance, in a Fermi glass the AC conductivity is expected 
to be a power law a oc ([55]! . In a Coulomb glass the 
conductivity is expected to take the form a (x U{r) ■u+uj'^ 
([S7)l. Here, the non-interacting Fermi glass functional 
form is returned at high frequencies, but a linear depen- 
dence is found at low frequencies. The crossover between 
these regimes is smooth and set by U{r) which is the typ- 
ical long range interaction strength between the charges 
forming the absorbing resonant pair. There are expected 
to be logarithmic corrections to these power laws. As the 
energy scale U{r) characterizes the long range Coulomb 
interaction, it is expected to go to zero at the MIT. As one 
approaches the quantum critical point, the behavior may 
also be interpreted using scaling laws (discussed below) 
that connect the temperature, frequency and concentra- 




FIG. 24 Absorption coefBcient in P doped Si normalized to 
the number of donors no for the three different labeled doping 
levels at 2 K. At the lowest doping, absorptions due to indi- 
vidual donor atoms can be seen. These absorption broaden 
into a impurity band at higher doping s. T he direct gap is at 
much higher energies as shown in Fig. 



25 From Ref. 
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FIG. 25 Optical conductivity of silicon doped close to the 
3D metal-insulator transition. The optical band gap is visible 
via the high energy onset at 1.12eV. The intermediate energy 
conductivity is associated with excitations from the dopant 
band to the conduction band. The conductivity at the lowest 
energies, i.e. photon assisted hopping conductivity, is due to 
intra-dopant band excitations. From Ref. H65|) . 



tion dependence of the response. It is expected that there 
is a correspondence between the frequency and tempera- 
ture dependent conductivity on both sides of the critical 
concentration. Such an analysis of the conductivity leads 
to a universal scaling function and defines critical expo- 
nents as discussed below. 

In Fig. 26 the pioneering data of Stutzmann and Lee 
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(|3ip on Si:B demonstrates the real part of the conduc- 
tivity measured with a broad band microwave Corbino 
spectrometer. The data shows a remarkable observation 
of conductivity linear in frequency at low frequencies and 
a sharp crossover at U (r) to oj'^ at high frequencies. Work 
on Si:P on samples farther from the metal-insulator tran- 
sition shows the same behavior ([Ml albeit the en- 
ergy scale of the crossover is found at higher energies as 
it is set by the long range Coulomb interactions (which 
increases away from the MIT). Although the proposed 
functional form of Efros and Shkovlkii approximately fits 
the data of Stutzmann and Lee ([5T| (and of Helgren et 
al. dSUl), it fails to account for the very sharp crossover 
between power laws. Stutzmann and Lee (|3T]l interpreted 
the sharp crossover as deriving from a sharp feature in 
the density of states - the Coulomb gap (a depression at 
the Fermi energy in the density of states that is a conse- 
quence of the long range Coulomb interaction). However, 
such an interpretation is at odds with theory and more- 
over, was not consistent with the work of Helgren et al. 
(f64| . where it was shown that the crossover was set by 
the typical interaction strength. The sharp transition 
between power laws remains unexplained; although one 
may consider that this sharp crossover may arise from 
collective effects in electron motion (69) . 



Si:B 




FIG. 26 The real part of the complex frequency dependent 
conductivity for two samples of Si:B at 85% and 88% of the 
way to metal-insulator transition (31). Data has been mea- 
sured up to 20 GHz with a Corbino geometry microwave 
probe. The crossover in the frequency dependent conduc- 
tivity from linear to quadratic is qualitatively consistent with 
a crossover from Coulomb glass to Fermi glass behavior. 



The theoretical functional form and experimental ob- 
servation is remarkable, because the ordering of regimes 
is at odds with the typical situation in solid state physics. 
Typically it is that non-interacting functional forms are 
expected at low frequency and temperature, while the 
explicit effects of interactions are seen at higher fre- 
quencies. The canonical example of this ordering is the 
Fermi liquid, where it is the low frequency spectra which 
can be typically modeled in terms of quasi-free non- 
interacting particles (perhaps with renormalized masses). 



The higher energy spectra are typically more complicated 
and show the non-trivial effects of electron-electron and 
electron-boson interactions. Similar physics is exhibited 
in the heavy-fermion compounds in that it is the low fre- 
quency regime that can be understood in terms of con- 
ventional Boltzmann transport of heavy electrons. In 
contrast, in electron glasses, it is the low frequency spec- 
tra which shows the effect of interactions while the high 
frequency spectra returns the non-interacting functional 
form. 

This was not the original expectation. Anderson orig- 
inally coined the term "Fermi glass" in analogy with 
the Fermi liquid, to describe a ensemble of localized 
charges whose properties where largely determined by 
Fermi statistics alone (fTOjl . 

" (The) Fermi liquid theorem is a rigorous 
consequence of the exclusion principle, it 
happens because the phase space available 
for real interactions decreases so rapidly 
(as or T^) . The theorem is equally true 
for the localized case: at sufficiently 
low temperatures or frequencies the 
non-interacting theory must be correct, 
even though the interactions are not 
particularly small or short range: thus the 
non-interacting theory is physically correct: 
the electrons can form a Fermi glass." 

- P.W. Anderson, 1970 

It is surprising that even in a material as thoroughly 
studied as doped bulk silicon that there exists no clear 
consensus as to the ground state and the nature of the 
low energy excitations at low dopings. It appears how- 
ever, that Anderson's speculation is not borne out by the 
experimental situation. The low lying excitations appear 
fundamentally changed by interactions. Of course, ex- 
isting experiments do not preclude that at even lower 
energies the Coulomb glass-like behavior will not break 
down or give way to a different response like a Fermi 
glass. 

D. Mott insulators 

Mott insulators differ from band insulators and Ander- 
son insulators in that their insulating effects derive from 
correlations and not from filled bands or disorder-driven 
localization. One expects to see very interesting aspects 
of spectral weight transfer from high energies to low as 
charges are doped into them and they eventually become 
metals. 

At half-filling it is believed that such Mott insulators 
are characterized by upper and lower Hubbard bands, 
which arc split by an energy U that is the energy to dou- 
bly occupy a single site. If one wants to occupy a site 
doubly, one must pay this energy cost U. As will be dis- 
cussed below the spectral weight J a{Lo)duj is a conserved 
quantity which depends on the total amount of charge in 
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a system. It is interesting to account for how spectral 
weight is transferee! from high to low energies by doping 
electrons into such a system ([7T|) . 



PES ^— E|, — IPES 



At half-filling, the lower Hubbard band is filled (oc- 
cupied) and the upper Hubbard band is unfilled. Since 
every site in the system has one electron, the spectral 
weight of upper and lower Hubbard bands are equal. Now 
consider the situation after one dopes a single electron 
to the upper Hubbard band. The unoccupied spectral 
weight of the upper Hubbard band reflects the total num- 
ber of sites which can have an electron added, the lower 
Hubbard band reflects the total number of sites occupied 
by a single electron, and the occupied weight of the upper 
Hubbard band reflects the total number of sites that are 
doubly occupied. These considerations then mean that 
doping a single electron transfers one state each from 
both upper and lower Hubbard bands as shown in Fig. |27] 
because doping changes the number of sites which are oc- 
cupied/unoccupied. It is natural to expect that a state is 
removed from the unoccupied part of the upper Hubbard 
band as now a state is occupied. However a state is also 
removed from the occupied part of the lower Hubbard 
band as its spectral weight quantifies how many elec- 
trons sit on singly occupied sites (and hence how many 
can be removed from singly occupied sites). The occupied 
weight of the upper Hubbard band is then '2', and the 
spectral weights of occupied lower Hubbard band and 
unoccupied upper Hubbard band are both N-1. The sit- 
uation differs from doping a band semiconductor because 
then doping a single electron into the upper band, only 
takes that single state from the upper band. There is no 
transfer of spectral weight across the gap. An interme- 
diate case is expected for the charge-transfer insulators, 
which the parent compounds of the high Tc supercon- 
ductors are believed to be. In these compounds a charge 
transfer band (consisting primarily of O 2p states) sits 
in the gap between upper and lower Hubbard bands and 
plays the role of an effective lower Hubbard band. The 
energy to transfer charge from from an oxygen to cop- 
per A becomes the effective onsite Hubbard U. Issues of 
spectral weight transfer in Hubbard and charge transfer 
insulators are discussed more completely in Ref. (j7ip . 



In Fig. [28] we show the doping dependence of the room 
temperature optical conductivity of La2_a;Sra.Cu04. At 
X = the spectra show a clear charge transfer gap of 
about 1.8 eV. This largely reflects excitations from O 
2p^ states to Cu 3d^^. Upon doping with holes, spectral 
weight is observed to move from the high energies to low 
in a manner consistent with the picture for the charge 
transfer/Mott Hubbard insulators in Fig. 



27 It is in- 



teresting to note that a remnant of the charge transfer 
band remains even for samples that have become super- 
conductors and are known to have a large Fermi surface. 
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FIG. 27 A schematic of the electron-removal and addition 
spectra for a simple semiconductor (left), a Mott-Hubbard 
system in the localized limit (middle) and a charge transfer 
system in the localized limit, (a) Undoped (half-filled) (b) 
one-electron doped, and (c) one-hole doped. The onsite re- 
pulsion U and the charge-transfer energy A are indicated. 
From Ref. (|7T|l . 
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FIG. 28 Doping dependence of the room temperature optical 
conductivity of La2-iSra:Cu04. At a; = the spectra show 
a clear charge transfer gap of about 1.8 eV. Upon doping 
with holes, spectral weight is observed to move from the high 
energies to low in a manner consistent with Fig. |27[ Figure 
adapted from (|72|) . 



E. Superconductors and other BCS-like states 

One can get a rough intuition of the electrodynamic 
response of superconductors for Huj ^ 2A and T^Tc 
from the relations in Eq. |25| which give the dissipation- 
less limit for the Drude model. One expects a S function 
peak at zero frequency in the real conductivity and a 
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l/w dependence in the imaginary part, whose coefficient 
is set by the strength of the S function. At frequencies on 
the order of the superconducting gap 2A, Cooper pairs 
can be broken and enhanced absorption should be found. 
At higher temperatures, thermally excited quasiparticles 
are created and one expects that as they are subject to 
essentially normal state dissipative processes, they will 
give a contribution to ai at finite lu. 

Very roughly one expects a two fluid scenario where 
the conductivity is approximately given by <t{uj) = 

|^(5(a; = 0)+i^^^^+ain{u!,T) where (Ti„ is a normal 
fluid component. At low temperatures, one expects that 
the superfluid density N{T) is degraded as 1 — e"'^/'^^^, 
as Cooper pairs thermally disassociate. As the superfluid 
density decreases the normal fluid component increases 
and one may in some circumstances expect an increase 
in ain to go as ^-^^I^bT ^ 

These above considerations are only approximate how- 
ever. The explicit temperature dependence of the gap 
and the mutual screening effects of the superfluid and 
normal fluid must be taken into account. Among other 
things, this means that the functional form of the contri- 
bution from 'normal' electrons cr„, although still peaked 
at w = 0, will be decidedly non-Lorentzian. Within the 
context of the BCS theory the response of a superconduc- 
tor can be calculated from the Mattis-Bardeen formalism. 
I will not go into the details of the calculation here. In- 
terested readers should consult Tinkham (5J . One of the 
important results from this formalism however is the ef- 
fect of the superconducting coherence factors, which de- 
pending on their sign can lead to enhanced or suppressed 
absorption over the normal state above and/or below the 
gap edge. This means that different symmetries of the 
superconducting order parameter (or order parameters 
of other BCS-like condensates like spin- or charge-density 
waves), leave their signatures on the dissipative response. 




FIG. 29 Temperature and frequency dependence of the con- 
ductivity o"! as evaluated from the Mattis-Bardeen expres- 
sion. The coherence peak exists only at low frequencies 
hjjj/2A < 0.1. The cusp in the surface corresponds to the 
energy gap. From Ref. (f75)) . 



Shown in Fig. 29 is the temperature and frequency 
dependence of the real part of the conductivity txi for 
uj > as evaluated from the Mattis-Bardeen expression 
for 'type IF coherence factors, which are appropriate for 
s-wave superconductivity. One can see that at T=0, one 
has no below gap absorption, and then a gentle rise be- 
gins starting at the gap edge 2 A. Despite the singularity 
in the density of states, the conductivity is suppressed at 
the gap edge, as the type II coherence factors effectively 
cancel the enhanced density of states. As one warms 
the sample, below gap excitations are possible from ther- 
mally excited quasiparticles. Just as above the gap edge 
the absorption is suppressed, the absorption is enhanced 
over the normal state at energies below the gap. This is 
again the consequence of type II coherence factors. 

This anomalous behavior is more clear along the tem- 
perature axis. If one follows the conductivity at a par- 
ticular frequency as a function of temperature, one will 
see an enhanced absorbence at temperatures below T^. 
This result although surprising is direct evidence for co- 
herence effects in the superconductor. An analogous en- 
hancement found in the nuclear relaxation rate via NMR 
- the so-called Hebel-Schlichter peak - was early evidence 
for the BCS theory (l73|) . 

These features contrast with the situation in Type I 
coherence factors as shown in Fig. |30) In a Type I co- 
herence factor SDW material one expects a very different 
response. There the coherence factors don't cancel the 
singularity in the density of states. Likewise, there is no 
Hebel-Schlicter-like peak in the temperature dependence. 
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FIG. 30 (Color) (left) Approximate frequency dependence of 
the dissipation rate from the Mattis-Bardeen theory for Type 
I and Type II coherence factors, (right) Approximate tem- 
perature dependence of the dissipation rate for u) — 0.1 A for 
both coherence factors. 



Qualitatively the behavior for Type II coherence fac- 
tors is realized in conventional s-wave superconductors. 
In Fig. |31| the normal state normalized conductivity 
from some of the original measurements of Palmer and 
Tinkham (|74p is shown, which is in good agreement with 
the Mattis-Bardeen prediction. The suppression at the 
gap edge is consistent with type II coherence factors and 
hence a s-wave superconducting state. Mattis-Bardeen 
predicts that subgap absorptions rise with increasing 
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temperature and become enhanced over the normal state 
conductivity. 
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FIG. 31 Frequency dependence of the normalized conductiv- 
ity through lead films at three different temperatures. The 
results are obtained through a combination of reflection and 
transmission measurements. The solid line is a calculation 
using the Mattis-Bardeen theory with a 2A/ft — 22.5 cm~^. 
Tc for these lead films is 7.2 ± 0.2 K. From Ref. (,74,) . 



In Fig. 



32 the conductivity of Nb at 60 GHz as mea- 
sured in a microwave resonance cavity is shown as a 
function of temperature (TZSjl . The peak below Tc is the 
electromagnetic equivalent of the Hebel-Schlicter peak in 
NMR. This is again consistent with s-wave superconduc- 
tivity. 

As noted, different behavior is expected for type I co- 
herence factors. There one expects an absorption en- 
hancehment at the gap edge and no 'Hebel-Schlicter' 
peak. The lack of such a peak in electromagnetic ab- 
sorption of cuprate superconductors is at least partial ev- 
idence for d-wave superconductivity in those compounds. 
It is due in part to the lack of a strong singularity in the 
d-wave density of states, but also because the coherence 
factor vanishes for g = tt, tt since A^, w — Aj._|.7r,7r for k 
near 7r,0. Such data is shown in Fig. |33| where the con- 
ductivity measured at a number of different frequencies 
in the GHz range on YBCO is displayed as a function of 
temperature (j76p . The data shows a broad peak in g\ at 
around approximately 35 K, but a comparison with the 
Mattis-Bardeen prediction shows that it is incapable of 
describing it as the rise in a\ is much more gradual than 
predicted. An explicit comparison is shown in Fig. |34[ 
This peak has been quite reasonably described in terms 
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FIG. 32 Temperature dependence of the complex conduc- 
tivity (t{lo) of Nb as evaluated from surface impedance mea- 
surements at 60 GHz in cavities. The solid curve is the weak 
coupling Mattis-Bardeen prediction. The dashed curve is a 
strong coupling Eliashberg prediction. From Ref. H75p . 



of a collapsing scattering rate below Tc of quasiparticles 
whose number remains appreciable until low tempera- 
tures due to the d-wave nature of this compound. An 
examination of the Drude model shows that in general 
one expects a peak in ui when the probe frequency is ap- 
proximately equal to 1/r the scattering rate. The exhib- 
ited peak is therefore not a Hebel-Schlicter peak, and the 
data is consistent with d-wave superconductivity. These 
experiments give evidence for d-wave superconductivity 
twice: the lack of a Hebel-Schlicter peak and the large 
CTi due to quasiparticle effects. Note that the other sharp 
peak right at T^ is believed to be due to fluctuations of 
superconductivity. Its functional form is also inconsistent 
with it being a Hebel-Schlicter peak. 
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FIG. 33 The real part of the conductivity extracted from 
the microwave surface resistance of YBCO at a number of 
different GHz frequencies. From ([76}. 



For completeness, I should mention that spin density 
wave compounds which can be treated within the BCS 
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FIG. 34 The real part of the conductivity extracted from the 
microwave surface resistance of YBCO. The BCS conductivity 
(sohd) is calculated using a Tc of 91.8 K, a gap ratio of 3.52, 
and various other physical parameters such as penetration 
depth, coherence length, and mean free path. It is completely 
incapable of describing the EM response, which shows the d- 
wave nature of these materials. From Ref. (|77|). 



formalism are expected to have a type I order parameter 
and a gap edge absorption enhancement. In many cases 
they can also be treated within a BCS formalism. See 
Ref. ()75)l for exphcit details. 

A vast literature exists on the electrodynamics of su- 
perconductors. I have given only the most superficial 
treatment here. Many considerations go into how super- 
conductivity is exhibited in optical spectra. See Tinkham 
([5), for instance, for in-depth discussions on the clean and 
dirty limits of superconductivity and Basov and Timusk 
(j9|) for the state-of-the-art on high-temperature cuprate 
superconductivity. 



V. ADVANCED ANALYSIS 
A. Sum Rules 

The optical constants satisfy a number of different sum 
rules (|3|). The most frequently used sum over the real 
conductivity is the one that will be discussed in detail 
here. As alluded to above, there is a relation for the 
integral over all frequencies of the real part of the con- 
ductivity to the total number of charges n and their mass 

Trie- 

f°° 47rne^ 
8 / (TiHdc^= (51) 

Jo "le 



Analysis of data in terms of sum rules provides a pow- 
erful tool that can be used to study spectral weight dis- 
tributions in a relatively model-free fashion. The above 
integral which extends from zero to infinity is the global 
oscillator strength sum rule, which relates the integral of 
the (Ji to the density of particles and their bare mass. It 
is related to the oscillator strength fs_s> discussed above 
in the context of the Kubo formula. Both are manifesta- 
tions of the /-sum rule of elementary quantum mechan- 
ics. Sum rules follows as consequence analyticity of the 
response functions and causaHty. See Refs. (HJ |3| for a 
few derivations. 
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FIG. 35 The effective number of carriers ne//(nc) as a func- 
tion of cutoff frequency Qc for Aluminum. A number of dif- 
ferent sum rules are investigated and displayed. The sum 
rule for €2 is equivalent to the sum rule expressed in Eq. |51| 
Figure adapted from Ref. (|79|l . 



In real systems the integral to infinite frequencies is 
rarely utilized in practice and we are usually more con- 
cerned with partial sum rules such as 

8 / cTi{oj)dLo= —. (52) 

Jo mb 

where W is the unrenormalized electronic bandwidth and 
rUb is the band mass. Sum rules of this fashion can be 
used to get information about occupation and correlation 
in a subset of bands. Of course if the integral cut-off W 
is extended to high enough energies the full sum rule Eq. 
|51| must be fufilled This puts a fundamental constraint 
on the various values that the individual sub-band band 
masses can assume. A partial sum rule analysis of the 
equivalent quantity 62 applied to the simple metal case 
of aluminum is shown in Fig. |35[ Aluminum's nominal 
electronic configuration is ls^2p^2p^3s^3p^ . Depending 
on the cutoff of the integral, the sum rule for different 
orbitals is satisfied. Integrating up to approximately the 
plasma frequency (15.5 eV) returns the number of elec- 
trons in the n = 3 valence band. At the x-ray L edge at 
around 100 eV, one begins to reveal spectral weight of 
the n = 2 formed bands. One can see that if the integral 
is performed to high enough energies (well above the K 
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edge of aluminum) one recovers the aluminum's atomic 
number '13' if one assumes the free-electron mass in Eq. 
[52] It is interesting to note that the rate of convergence 
is different for the different sum rules plotted. This is re- 
lated to how spectral is distributed in the spectra and can 
be related to quantities like the strength of interactions. 

In correlated systems one frequently makes use of an 
even more limited sum rule and only performs the inte- 
gral over an energy only several times the Drude width 
to determine the mass renormalized due to interactions. 
For instance the small spectral weight in the very nar- 
row Drude peak as shown in Fig. [8] gives a very large 
mass (~ lO^mg). An integral as such over a small renor- 
malized Drude peak with small spectral weight is how 
the very large masses in Fig. |4] were generated. One 
can view this small mass as coming from a very narrow 
interaction derived band where the integral only has to 
be performed over a small fraction of W to capture its 



spectral weight. However, the sum rules of Eqs. 51 and 
|52]must ultimately be satisfied at high enough energies. 
This means that there must be be some higher energy 
satellite in ai that contains enough spectral weight to 
return the band mass rrib when the integral in Eq. |52] 
is performed out to the unrenormalized band width W. 
Note that similar information can be gained from apply- 
ing the extended Drude model analysis, which is outlined 
in the next section (Sec. V.B). 

Sum rules have been used extensively in the analysis of 
the data for correlated electron systems because, in some 
circumstances, they allow a relatively model independent 
way to analyze the data. For instance, one can show that 
for a tight binding model with nearest neighbor hopping 
only that 
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(Ti (w) duj 



(53) 



where is the lattice Constance in the incident E field's 
polarization direction and Kr is the effective kinetic en- 
ergy. This relation has been used extensively in the 
cuprate superconductors to attempt to show evidence for 
or against a novel lowering of the kinetic energy-driven 
mechanism for superconductivity. This is in contrast to 
the usual mechanism in the BCS theory where it is the 
potential energy that drops when falling into the super- 
conducting state. Such analysis is not trivial in these 
materials, as there is evidence that the spectral weight 
that goes into the condensate is derived not just from the 
near Ep electronic band over a width much smaller than 
W , but over an energy range many times the band width 
W 
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FIG. 36 Temperature dependence of the low-frequency spec- 
tral weight Ai+_d(T) and the high-frequency spectral weight 
Ah{T), for optimally doped and underdoped Bi2212. Insets: 
derivatives -T'^dAi+o / dT and T'^dAih/dT. From Ref. 



It is believed that the usual sum rule for superconduc- 
tors (Ferrell-Glover-Tinkham) (5) is satisfied for the ab 
plane conductivity of the cuprate superconductors up to 
about 10% accuracy (9). Molegraaf et al. ^80) used a 
combination of reflectivity and spectroscopic ellipsome- 
try in Bi2212 (Fig. [36| to conclude that about 0.2 - 0.3% 
of the total strength of the superconducting i5 funtion 
is collected from an energy range beyond 10,000 cm"^ 
to 20,000 cm~i and that the kinetic energy decreases 
in the superconducting state. This is a very small effect, 
but still large enough to account for the condensation en- 
ergy. Such analysis of total spectral weight transfer can 
be tricky as it requires the analysis of spectral weight be- 
yond the measurement regime, which is only constrained 
through the Kramers-Kronning relations. These results 
or rather the analysis of the data has been disputed by 
Boris et al., (|8ip who form a nearly opposite conclusion 
regarding spectral weight transfer based on nearly iden- 
tical experimental data. 

The c-axis response has been investigated in a simi- 
lar manner and a similar transfer of spectral weight from 
very high energies to the S function has been observed, 
although in this case the magnitude of the effect is defini- 



It may be that this extreme mixing of low and high energy scales 
between, which prohibits a true low energy description is an al- 
ternative definition of a correlated material. 



10,000 cm~^ is approximately equivalent to 1.24 eV and 14,390 
K! It is remarkable that such energy scales are explicitly relevant 
for a 100 K phenomenon. In contrast conventional superconduc- 
tors satisfy the 'Ferrell-Glover-Tinkham' sum rule at energies 
approximately 10 - 20 A (O. 
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tively insufficient to account for Tc 

As mentioned above, there are a number of sum rules 
in addition to the most commonly applied one for cti 
(or equivalently the first moment sum of £2)- Prob- 
ably next most commonly applied is the sum rule for 
the first moment of the loss function, whose integral is 
equal to that for the first moment sum of £2 ®. This 
has been called the longitudinal /-sum rule. A simi- 
lar sum rule concerns the first moment sum of the ab- 
sorption coefficient k{uj) (79]). There is also a sum rule 
on the loss function itself, that can be related to the 
Coulomb correlation energy of the ground state wave- 
function (f5^ Similarly, there is second moment 
conductivity sum rule and a third moment loss func- 
tion sum rule that in certain circumstances allows one 
to evaluate the strength of Umklapp scattering and the 
static electron-lattice energy ((Ml [S5J. There are also 
first inverse moment sum rules for the conductivity and 
loss function, which are equal to various combinations 
of pure numbers (For instance ~-^™ ^£(a)) '^^ ~ 
These allow the calibration of, for example, the absolute 
intensity of an electron energy loss measurement since 
these sums have no dependence on any material param- 
eters. Other related sum rules include the inertial sum 
rule on the index of refraction (J^ [n{uj) ~ l]duj ~ 0) 
(IHSl [53 [HHI and the "DC sum rule" on the dielectric 
function {J^[€i{uj) - l]duj = -2n^aDc) dHUHZl)- There 
is also a sum rule on the frequency dependent scattering 
rate ()95| . which is discussed below. For further details 
on sum rules, see the very extensive literature on this 

subject (I3l3|7ilHllMl[Hl[H51[HSl[HZllHilMl)- 



B. Extended Drude Model 



Within the Drude-Lorentz (Sec. II. B ) model we con 



sider that conduction electrons are quasi-free. The scat- 
tering rate of the Drude intraband contribution is con- 
sidered to be characterized by a single frequency inde- 
pendent relaxation time r. It is well known however that 
various inelastic channels can contribute to relaxation in 
solids, each with a characteristic frequency dependence. 
For instance, in 3D the fermion-fermion scattering rate is 
expected to exhibit an dependence at low frequency. 
An electron scattering off an Einstein boson (a single 
spectral mode with a well-defined frequency luq) is ex- 
pected to show an onset in scattering at the boson fre- 
quency cjo which is the threshold to make a real (as op- 
posed to virtual) excitation. 

Allen and Mikkelsen (.89) proposed that one could 
capture these frequency dependences in optical spectra 
through an extended Drude model (EDM) where the 
mass and scattering rate are explicitly frequency depen- 
dent. Inverting the complex conductivity one gets 
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where rrii, is the band mass. 

Both in spirit and formalism this is equivalent to 
standard treatments in many-body physics where one 
posits that the effects of interactions can be captured 
by shifting the energy of an added electron from the 
bare non-interacting energy e(k) by a complex self- 
energy E(a;,k) = Ei(ti;,k) + iTt2{uj ,'k.) ■ One can express 
the optical response in terms of the complex self-energy 



1^2^(0;) as 
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and in terms of previously defined quantities 
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A caveat should be given about the interpretation of 
optical self-energies. Although an identification can be 
made with the quasi-particle self-energies that appear 
in one-particle spectral functions, the two self-energies 
are not exactly the same. While quasi-particle scatter- 
ing depends only on the total charge scattering rate to 
all final states, transport and optics experiments have a 
much larger contribution from backward scattering than 
forward scattering as it degrades the momentum much 
more efficiently. In general this means that the quan- 
tities derived via optics emphasize backward scattering 
over forward scattering and so the optically derived quan- 
tities can contain a vertex correction not included in the 
full quasi-particle interaction. In this sense the two self- 
energies are formally different although they can contain 
much of the same information. 

Frequently the extended Drude model will be in- 
troduced in terms of a complex frequency dependent 
Memory function, M^uj) = l/r(aj) — iLuX^u), where t 
is the lifetime and 1 -I- A = m* /mb 
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One can put Eq. [59] in the form of the standard Drude 
model with the substitutions t*(w) = [1 -I- A(aj)]r(a;) and 
= [1 + to get 
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This equation describes the frequency dependence of a 
particle with renormalized plasma frequency (and hence 
renormalized mass) and renormalized scattering rate. 
Compare Eq. |60] with Eq. |26) Note that the quantity 



1/r* is not equivalent to the quantity defined in Eq. 55 
above, as 1 /t* includes the renormalization effects of the 
lifetime as well as the mass, whereas 1/r includes only 
lifetime effects. The two quantities differ by a factor of 
1 + A(w). In this sense one may argue that 1/t{lo) is the 
more intrinsic quantity, as also evidence by its direct pro- 
portionality to the imaginary part of the self-energy, but 
this point is essentially a semantic one. l/r* is the effec- 
tive scattering rate (or dressed scattering rate) however, 
because for weak frequency dependence it is the actual 
half width of the Drude peak in the optical conductivity. 

Various contributions to the total scattering rate such 
as electron-electron (ee), electron-phonon (ep), electron 
magnon (em), etc. add within the prescription set by 
Matthiesen's rule as ^ = — + — + — + Note 

t{u)) To Te,, Tcp Tern 

that there is typically also a contribution to the elas- 
tic scattering which is frequency independent that comes 
from static disorder. As mentioned above this can be 
seen as a measure of the degree to which the translational 
symmetry is broken by disorder and hence the degree to 
which strict momentum conservation in the optical selec- 
tion rule q — Q can be violated. 

As -^r^j^j^ and m*{ijj) are a particular parameterization 
of the real and imaginary parts of a complex response 
function, they are Kramers-Kronig related. Analogous 
to the case of a(uj) detailed above, if one knows —t-^ for 
all uj, then m*(u}) can be calculated and vice versa. 

The near-Ei? electrons of the noble metals, like Cu, 
Ag, and Au have free-electron-like s-states, which form 
a Fermi surface that deviates from almost spherical sym- 
metry only near the Brillouin zone necks in the < 111 > 
directions. The above analysis applied to such metals 
is shown in Fig. |37] and Fig. |38j For Cu and Ag the 
scattering rate shows an approximately dependence of 
the scattering rate, which is in accord with the expecta- 
tion for electron-electron scattering. It is interesting that 
electron-electron scattering is observed up to such high 
energies as we may have expected other contributions 
such as electron-phonon scattering also to give an ob- 
servable contribution. We should reiterate that a purely 
translational invariant electron gas, cannot dissipate mo- 
mentum by electron-electron collisions alone. In this re- 
gard effects like umklapp or interband (Baber) scattering 
are essential to see the effects of such collisions in the op- 
tical response. Also note here the difference between the 
bare scattering rate 1/r and the dressed scattering rate 
1 /r* . The difference is only significant when there is a 
large renormalization in the effective mass. 

Note that in Eqs. [54| and [55] there is the matter of how 
to define the plasma frequency LOp . For the electron 



gas the plasma frequency is given by Aivne^ /m, however 
the precise definition is less clear in a real material with 
interactions, interband transitions, and orbital mixing. 
In the above analysis must come from the spectral 
weight of the full intra-band contribution. Rigourously 
the plasma frequency should be defined from the sum 
rule Eq. [STjas 



(61) 



where contribution to the conductivity com- 

ing from all (and only) intraband spectral weight . 
The integral extends all the way to infinity, which gives 
of course a practical difficulty in real systems because 
the value of the spectral weight will be contaminated by 
interband contributions. Note that for a complicated in- 
teracting system the integral and ujp must include not 
only the contributions of the weight of the low frequency 
Drude peak, but must also include higher energy parts 
of the spectra (satellites) if they derive from the same 
non-interacting band. It can be difficult experimentally 
to estimate ujp accurately. For instance, in the case of 
a system with strong electron-phonon coupling the spec- 
tra may be approximately modeled as the usual Drude 
model with a simple co — Lorenztian and frequency 
independent scattering rate at low frequencies, but with 
a higher frequency satellite at the characteristic phonon 
frequency which results physically from the excitation of 
a real phonon as well as an electron-hole pair. The high 
frequency satellite manifests in the scattering rate as an 
onset at the phonon frequency as an additional decay 
channel becomes available at that energy. In this sense 
the plasma frequency which goes into an extended Drude 
analysis is the plasma frequency associated with the spec- 
tral weight of both the low frequency 'Drude' contribu- 
tion and the satellite. The spectral weight of the low 
frequency part may be clear, but it may be a difficult 
practical matter to clearly identify the spectral weight of 
the satellite, due to the overlapping contributions of it 
with true interband transitions etc. In practice, a high 
frequency cutoff to the integral in Eq. [61] is usually set 
at some value which is believed to capture most of the 
intraband contribution, while minimizing contamination 
by interband transitions. A number of different criteria 
can be used: the requirement that the integral is temper- 
ature independent, an appeal to reliable band structure 
calculations etc. etc. 



Note that this is only important for the overall scale factor of 



these expressions. Their functional dependence and form is in- 
dependent of the choice of the plasma frequency. 
Note that in anything other than a simple Drude system the 
plasma frequency here is not given by the plasma minimum in 
the reflectivity as that quantity is renormalized by a factor of 
^1 — As mentioned above, this later quantity is known as the 
screened plasma frequency or plasmon frequency. 
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FIG. 37 Scattering rate vs. oj^ of (a) Cu, (b) Ag and (c) Au. Solid and dotted lines are the bare scattering rate l/r and the 
dressed scattering rate l/r* respectively as detailed in the text. From (|90|l 
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FIG. 38 Mass enhancement factor 1 + A of noble metals. From 



There can also be a difficulty in definitions of the 
plasma frequency for systems described by the Mott- 
Hubbard model, where a strong onsite repulsion splits 
a single metallic band into upper and lower Hubbard 
'bands'. In such cases, one has one band or two depend- 
ing on definitions and so the matter of the true intraband 
spectral weight can be poorly defined. In such cases it is 
important to remember that in this and other cases that 
the Drude model is a classical model applicable only to 
weakly interacting mobile particles. Although the EDM 
can be put on more rigorous ground expressing its quan- 
tities as optical self-energies, one still is always making 
an analogy and connection to a non-interacting system. 
To the extent that the EDM is just a parameterization 
of optical spectra, it is always valid. However, to the 
extent that its output can be interpreted as a real mass 
and scattering rate of something (i.e. a particle), it is 
important that that something exists. In other words, it 
must be valid to discuss the existence of well-defined elec- 
tronic excitations in the Fermi-liquid sense in the energy 
regimes of interest. 

In this regard, one must also be aware of the contri- 
bution of any interband terms in the evaluation of Eqs. 



for excitations that are nominally w = 0. In principle 
finite frequency interband excitations (classically Drude- 
Lorentz) can overlap with a broad Drude intraband term 
giving artifacts that should not be considered as contri- 
butions to the scattering rate and effective mass. The 



sharp onset in 1/t(w) in Fig. 37 and the cusp in Fig 
are due to a onset of strong interband excitations in t 
spectral region (Fig 
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The EDM has been used in many different contexts 
in strongly correlated systems. An ideal example is af- 
forded by its application to heavy fermion systems. As 
mentioned above, due to interaction of conduction elec- 
trons with localized moments the conduction electrons 
can undergo extremely large mass renormalizations of 
a factor of a few hundred over the free electron mass. 
These mass renormalizations develop below a coherence 
temperature T* . Such effects are reflected in the optical 
spectra. As seen in Fig. |39] at high temperatures an 
EDM analysis of CeAla reveals a pure Drude-like, essen- 
tially frequency independent mass and scattering rate. 
The conventionally sized mass and large scattering rate 
at high temperatures reflect the effects of essentially nor- 
mal electrons scattering strongly and incoherently with 
the localized moments The spectra show no partic- 
ularly interesting frequency dependence. Below the co- 
herence temperature at low frequencies significant renor- 
malizations to the mass and scattering rate are seen. The 
mass develops a significant enhancement and the scatter- 
ing rate drops. The coherence temperature T* and the 
frequency above which the unrenormalized values are re- 
covered ujc are equivalent to each other to within factors 
of order unity (with H = ks = 1). Remarkably the opti- 
cal mass as compared to the band mass is increased by a 



54 and 55 These expressions for the EDM are only valid 



Note that the conventional EDM analysis expresses m* as a ratio 
with the band mass mi, if u!p is treated correctly as discussed 
above. This is in contrast to other techniques like, for instance, 
de Haas- van Alphcn that expresses m* as a ratio with the free 
electron mass m». 
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factor of almost 350. Note that although both the mass 
and the low frequency scattering rate undergo dramatic 
changes, the DC Drude conductivity is not changed sig- 
nifiantly, as udc r/m and the respective changes ap- 
proximately cancel. This is reflected in the fact that the 
relation m*{uj 0)/mf, is expected to be equivalent to 
t/t*{lli 0) to within factors of order unity ([TT)| . Note 
that these changes in scattering rate and mass are re- 
flected in the optical conductivity itself by the formation 
of a very sharp zero frequency mode (essentially a Drude 
peak) which rides on top of the normal Drude peak (|lip . 

Evidence that the renormalizations in the optical mass 
reflect the existence of real heavy particles in these sys- 
tems and not just a convenient parametrization of the 
optical spectral can be seen in the above discussed Fig. 
[4] which shows a proportionality of the optically mea- 
sured mass to the linear coefficient of the specific heat 
(which is proportional to inverse of the near E^^ density 
of states and hence the mass). It is remarkable that a pa- 
rameter determined dynamically - essentially by shaking 
charge with an oscillating E-field - is precisely related 
to a quantity which is determined thermodynamically 
by quantifyng the amount of heat absorbed. This graph 
provides a remarkable demonstration of the quasiparticle 
concept! 
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FIG. 39 Scattering rate and effective mass in CeAla. (a) 
Frequency-dependent optical scattering rate for CeAIs at four 
temperatures, compiled using both infrared reflectivity and 
microwave cavity data, (b) Frequency dependence of m* ob- 
tained from Eq. [551 (193)) . 




veal themselves in the optical spectra. Spectral features 
in the optical conductivity deriving from electron-boson 
coupling were orginally discussed in the context of Hol- 
stein processes - the creation of a real phonon and an 
electron-hole pair with the absorption of a photon - in 
superconducting Pb films (jU [92}. The EDM analysis 
can reveal significant information about the interaction 
of electrons with various bosonic modes. 
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FIG. 40 Electron-boson model calculations with a boson 
spectral density A{Q,) = ujoS{u! — uo)- The top panel gives 
the optical conductivity and the lower two panels show the 
corresponding scattering rate and mass enhancement. The 
coupling constant in this calculation is set to 1. From l|94p . 



In Fig. |40] the optical conductivity calculated for an 
electron with impurity scattering coupled to a single Ein- 
stein phonon is shown for a few different temperatures. 
(See Puchkov, Allen, Shulga (94) for details regarding 
the calculation of the effects of electron-boson scatter- 
ing in optical spectral.) The lower two panels show the 
scattering rate and mass enhancement generated from 
the optical conductivity. At low temperatures the opti- 
cal conductivity (Fig. 40 (top)) reveals a sharp Drude 

which results from the usual 



peak at low frequencies, 
electron- hole excitations. The sharp onset in absorption 
at the phonon frequency is the threshold for excitation 
of a real phonon. At threshold, electron-hole pairs are 
excited additionally, but it is the phonon which carries 
away most of the energy. This is reflected in the scat- 
tering rate, which is flat and small at low to, but shows 
a cusp and sudden increase at loq as the phase space for 



As mentioned, in addition to electron-electron scat- 
tering, the effects of electron-boson interactions also re- 



scattering increases discontiously (Fig. 40 (middle)). 

Although real phonons can only be excited above 
threshold, electron-phonon coupling manifests itself at 
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energies below threshold by an increase in the effective 
mass (Fig. 40 (bottom)). Physically, this can be un- 



derstood as a renormalization of the electron-hole pair 
energy modified by the excitation of virtual phonons. 
This renormalization of the energy manifests itself as an 
increased mass. The electron-hole pair can be viewed as 
dragging a cloud of virtual phonons with it. 

One can get an approximate measure of the boson cou- 
pling function as 



2-K duj'^ ayu)) 



(62) 



Although the results of the calculation diplayed in Fig. 
[40] is for phonons, the general idea holds for electron- 
boson coupling generically (|95l [M]! although, of course, 
the bosonic density of states is generally more compli- 
cated than a single sharp mode. Recently Tediosi et al. 
(|97p found a remarkable example of this concept in semi- 
metal bismuth, where the bosonic excitation is formed 
out of the collective excitations of the electronic ensem- 
ble itself i.e. the plasmon. This is an effect that should 
be present in all metals, but is enhanced in a semi-metal 
like bismuth due to its very low carrier density and high 
Coo ~ 100. At low frequencies Tediosi et al. ()97) found 
in their EDM analysis (Fig. 41 ) that the mass and scat- 



tering rate are approximately frequency independent as 
expected for the Drude model. However, at higher fre- 
quencies there is sharp onset in the scattering rate at a 
strongly temperature dependent position. As one cools 
bismuth, the charge density changes dramatically and the 
plasmon frequency (as given by the zero crossing of ei) 
drops by a factor of almost two. It was found that the 
energy scale of this sharp onset in scattering tracks the 
independently measured plasmon frequency as shown in 



Fig. 42 One can note the strong resemblance of the data 
to the calculation of an electron interacting with sharp 
mode in Fig. |40] A strong coupling as such between elec- 
trons and plasmonic electron degrees of freedom may in 
fact be captured within the same Holstein Hamiltonian 
that is used to treat the electron - longitudinal phonon 
coupling and describe polarons. For this reason this col- 
lective composite excitation has been called a plasmaron 
([98] l99|l when exhibited in the single particle spectral 
function. The observation of Tediosi et al. was the first 
such observation optically. 

EDM analysis has been applied extensively to the high- 
Tc cuprate superconductors dS}. In Fig. 43 I show the 
imaginary and real parts of the optical self-energy as de- 
fined in Eq. [57] above for a series of four doping levels of 
the compound Bi2Sr2CaCu208+5 (|100p . The displayed 
samples span the doping range from the underdoped to 
the severely overdoped. Again one can notice the resem- 
blance of the data to that of the model calculation in 
Fig. 40 and indeed data of this kind, along with simi- 



lar signatures in angle-resolved photoemission have been 
interpreted as generic signatures of electron-boson inter- 
action in the cuprates. 
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FIG. 41 (Color) Frequency dependent scattering rate (a) and 
effective mass (b), calculated from the optical data of ele- 
mental bismuth using eoo. The low frequency scattering rate 
T~^{uj) progressively falls as the temperature is lowered. An 
approximately frequency independent region is interrupted by 
a sharp onset in scattering and a decrease of the effective 
mass, which was associated with plasmon scattering. 
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FIG. 42 (Color) A parametric plot of the sharp onset in 
scattering plotted versus the independently measured plas- 
mon frequency in elemental bismuth. The plasmon frequency 
changes monotonically as the sample is cooled from room tem- 
perature to 15K. 



The identity of this boson is a matter of intense current 
debate in the community. Some groups have pointed to 
a many-body electronic source (|96l llOOj llOlj I102p related 
to the '41 meV magnetic resonance mode discovered via 
neutron scattering (|103l I104p . Others have argued that 
these feature's presence above Tc, the relative doping in- 
dependence of its energy scale, and its universality among 
material classes demonstrate a phononic origin and in- 
dicates the strong role that lattice effects have on the 
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FIG. 43 (Color) a-d, The frequency and temperature depen- 
dent optical scattering rate, 1/t{uj) for four doping levels of 
Bi2Sr2CaCu2 08+i. a, = 67 K (underdoped); b, 96 K (op- 
timal); c, 82 K (overdoped); d, 60 K (overdoped). e-h, The 
real part of the optical self-energy, - 2E°''(aj) as defined in Eq. 
[57| From Ref. (fTDO)) . 



low-energy phvsics ()105p We note that whatever the 
origin of this effect, it is almost certainly the same as 
that which causes the upturn in the optical scattering 
rate d^l [Ml HOOP - We should also mention that much of 
the discussion regarding these self-energies and the like 
has implicitly assumed that the only possibility of under- 
standing such features are a consequence of some sort of 
electron-boson coupling, but this is not necessarily the 
case. There are proposals that purely electronic inter- 
actions and proximity to a Mott insulator can give such 



My (NPA's) opinion is that these are by and large straw-man 
arguments and perhaps not the right way to approach the prob- 
lem of the origin of high-Tc . In a strongly correlated system like 
the cuprates, it is clear that everything is strongly coupled to ev- 
erything else. Magnetism is strongly connected to lattice effects 
and vice versa. In this regard, it is not surprising that experi- 
ments sensitive to the lattice show lattice effects and experiments 
sensitive to magnetism show magnetic effects. One should not 
expect in such materials that various degrees of freedom neatly 
partition in separate distinct subsystems. 



features ([TTOl [TUTll . 



C. Frequency dependent scaling near quantum critical 
points. 

A quantum phase transition (QPT) is a zero- 
temperature transition between two distinct ground 
states as a function of a non-thermal parameter, such as 
magnetic field, pressure, charge density etc. (I108j I109p . 
Unlike conventional phase transitions, which are driven 
by thermal fluctuations, they are driven essentially by 
diverging fluctuations of the system's zero point motion. 
Once thought to be only of academic interest, the exis- 
tence of a QPT nearby in a material's parameter space is 
believed to influence a whole host of finite temperature 
properties. QPTs may hold the key to understanding 
the unusual behavior of many systems at the forefront of 
condensed matter physics pOSp . 

Just as in the case of conventional phase transitions, 
QPTs are believed to be characterized by diverging 
length and time scales if second order. In the disordered 
state, one envisions an order parameter that fluctuates 
slower and slower over longer and longer length scales as 
the transition is approached until at the transition these 
length and time scales diverge. These length and time 
scales are the correlation length and time Tc- The 
divergence of the fluctuation time scales, which in the 
case of conventional transitions is called critical slowing 
down implies the existence of a characteristic frequency 
ujc which vanishes at the transition. 

In a QPT, the transition occurs as a function of some 
non-thermal parameter K. It is usually considered that 
close to the transition at K^, diverges in a manner 
oc \K — Kc\~'^ and Tc diverges cx \K — Kd^^^" ^ where 
v and z are called the correlation length and dynamic 
exponent respectively. 

From the classical case a well established formalism 
exists for the understanding of various physical quanti- 
ties close to continuous (2nd order) transitions. Widom's 
scaling hypothesis assumes that close to the transition 
the only relevant length and time scales are those asso- 
ciated with diverging correlations in the order parame- 
ter piop . Such considerations hold for QPTs also and 
means that physical observables such as magnetization or 
conductivity are expected to have scaling forms, in which 
independent variables such as the temperature, probe fre- 
quency or wave vector appear in the argument only as a 
product with quantities like and Tc- This is expressed 
formally in terms of finite — size scaling. For an incident 
wave vector q — such a scaling form can be written 



OiK,Lu,T) 



:F{to/T,Tn). 



(63) 



Here O is some observable and F is the scaling func- 
tion, do is the scaling dimension, which determines how 
various physical quantities change under a renormaliza- 



34 



tion group transformation. It is often close to the "en- 
gineering" dimension of the system (for instance for 
2D conductivity), but can acquire an anomalous dimen- 
sion in certain circumstances. Generally one expects fre- 
quency and temperature dominated regimes of these re- 
sponse functions for uj ^ T,lOc and T 3> lo^lJc respec- 
tively. For further discussion on scaling forms and moti- 
vation on this and other related points see the excellent 
review by Sondhi et al. p08p . 

Scaling functions written as a function of and Tc 
can be very powerful as they allow the analysis of ex- 
perimental data independent of microscopic theory. For 
instance, the scaling exponents v and z only depend on 
certain global properties of the system such as symmetry, 
dimensionality, and the nature of the dominant interac- 
tions. 

There are many studies where the existence of a QPT 
is inferred from frequency dependent scaling (see for ex- 
ample (imp ) in a regime where frequency effects domi- 
nate over temperature ones, but there are far fewer stud- 
ies where scaling is investigated as one passes through 
a known QPT. I know of only two where scaling in the 
frequency dependent conductivity has been investigated. 
This is unfortunate because, as mentioned, such stud- 
ies could potentially be tremendously powerful for the 
investigation of correlated systems where the appropri- 
ate underlying physical model is not clear. This shows 
the extreme experimental challenges inherent in accessing 
the experimentally relevant frequency and temperature 
ranges where huj ^ ksT while both hio and ksT are as 
low as possible, but also below any higher energy scales 
that are not relevant to the ordered state. Hopefully this 
situation will change with the increasing prominence and 
feasiility of broadband microwave and THz techniques. 

The first experimental study to investigate frequency 
dependent scaling near a QPT was that of Engel et al. 
()112p , who used a waveguide coplanar transmission setup 
to perform broadband microwave spectroscopy on quan- 
tum Hall systems. In quantum Hall systems the lon- 
gitudinal resistivity goes to zero at field values where 
the transverse resistivity assumes perfectly quantized val- 
ues. Although the value of the Hall resistivity is perfectly 
quantized, the width of the transition region depends on 
temperature (|113p and on measurement frequency (|112p . 
In Fig. |44]the Re a^x vs. B is shown at three different 
frequencies and two different temperatures. One can see 
that at 50 mK, the width of the transitions regions (AS) 
is strongly frequency dependent, while at 470 mK the 
widths are almost frequency independent. 

This frequency and temperature dependence of AB 
(as defined as the extrema of the derivatives of Re a^x) 
is summarized in Fig. |45]for one of the spin-split levels. 
Similar data are found for other levels. One sees that for 
all such levels at low frequencies, the data are frequency 
independent and temperature dependent, whereas for low 
temperatures the data are temperature independent and 
frequency dependent. The crossover between regimes 
takes place when Shw ~ fc^T. Such data are consis- 
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FIG. 44 Re a^x vs. B at three frequencies and two tempera- 
tures. Peaks are marked with Landau level A*' and spin index. 
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FIG. 45 Peak width between extremal points in longitundinal 
resistance d Re Oxx/dB AB vs. measurement frequency at 
three different temperatures. Data are shown for the N = I 
Landau level. 



tent with scaling theories of quantum criticality, where 
the only relevant frequency scale at the QPT is set by the 
temperature itself. It is also significant that the crossover 
condition between frequency and temperature dominated 
regimes involves h, showing that the physics is essentially 
quantum mechanical. In the high frequency data it is 
observed that spin split levels give a, AB which is ap- 
proximately (X Lj'^ , with 7 = 0.43. This is consistent 
with scaling theories that give an exponent for quantum 
percolation of I/7 = ziy = 7/3 and also consistent with 
the temperature dependence of these transition widths 



35 



as measured by DC transport that show a depen- 
dence. It was claimed that through an analysis of both 
the temperature and frequency dependent exponents, one 
can extract the dynamic exponent z, which yields 2 = 1. 
This is consistent with long range Coulomb interaction. 

Both the DC data (|113p and the high frequency data 
have transition widths which exhibit a power law depen- 
dence on temperature or frequency with an exponent of 
3/7 at all spin-split transitions, showing that all such 
transitions between quantum Hall levels fall into the same 
universality class. It should be possible in such a study to 
incorporate both temperature and frequency dependence 
into a single scaling function although this has not been 
done yet. It is expected that a general scaling function 
of the form 



p{B,T,uj)^F{u/T,6/T^') 



(64) 



applies where 5 =\B — B^/Bc measures the distance to 
the quantum critical point. This function is equivalent 



to Eq. 63 using the fact that the scaling dimension of the 
resistivity vanishes in d = 2. In the high frequency and 
high temperature limits, Eq. |64| rcduccs to one where fre- 
quency or temperature appear in the argument as —n-^ 



and 



respectively. 
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FIG. 46 Re(7(aj) plotted against y/uj for temperatures 2.8 to 
32K and frequencies 87-1040 GHz. For the lowest temper- 
atures and when hu) > ksT, the data follow a y/uj. The 
data approaches the DC values plotted on the vertical axis 
for hu < kaT. 



The other frequency dependent study was that of 
Lee et al. (1114) I115P who applied the ideas of finite 
size scaling to the 3D metal-insulator (MI) transition in 
Nbi_a;Sia;. They measured the frequency dependent con- 
ductivity via millimeter wave transmission through thin 
films whose relative concentrations of Nb to Si were tuned 
exactly to the MI QPT point at Xc- As in the study of 
Engel et al. fll2() . for samples tuned right to quantum 
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FIG. 47 Log-log plot of conductivity scaled by the factor 



1/2N 



uj/T. For temperatures 16 K and 



below, the data for the entire experimental frequency range 
collapse onto a single curve within the experimental noise. 
Higher temperature data begin to rise above the collapsed 
data systematically for low scaled frequencies. The dashed 
line is a trial scaling function as described in the text. 



criticality it is observed that for the Hlu ^ ksT data 
the conductivity is temperature independent and that 
for Huj <C fcsT, the data are frequency independent (Fig. 
46 1 . A crossover from a frequency dominated regime to 
a temperature dominated regime implies the existence of 
a scaling function ^ of a form that can be inferred from 
above Eq. |63] 



(T(Xc,r, W) = CT' 



(65) 



In previous work on similar samples it was found that 
for samples tuned to the critical concentration the DC 
conductivity followed a T^/-^ relation implying z — 2. 



As shown in Fig. 47 successful scaling can be found by 



plotting the scaled conductivity Re a{T,uj)/CT'-'^ us- 
ing this z. This procedure collapses the data over the 
entire measured frequency range for temperatures 16 K 
and below, implying the applicability of a scaling func- 
tion that depends only on the scaled frequency. Data 
for temperatures higher than 16 K starts to rise above 
the other collapsed curves, indicating the appearance of 
other mechanisms (electron-phonon scattering for exam- 
ple) which limits the size of quantum fluctuations. This 
is consistent with the DC result where the conductivity 
starts to deviate from the T^ above 16 K. 
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As already inferred from Fig. [46] the scaling func- 
tion shows both temperature and frequency dominated 
regimes, with a crossover at approximately huj ~ kgT. 
Again, this shows that at the QPT the only energy scale 
is set by the temperature itself. Note, that this crossover 
happens at a slightly different oj/T ratio than the quan- 
tum Hall case, but is still of order unity as expected. For 
high frequencies the scaling function follows a power law 
dependence with the same exponent as was used to scale 
the vertical axis with z = 2. This means that an equiv- 
alent scaling function could have been found by divid- 
ing the conductivity data by a; 5 and plotting the data 
as a function of T/ui. Lee et al. attempted to 

guess the form of the scaling function with the expression 
^ = Re{\ — ibhu /ksTY/'^). The function successfully 
captures the high and low w/T, but misses the sharp 
crossover near hLu/ksT ^ 1. 

The value of the dynamic exponent z ~ 2, which was 
found in this study is interesting. It is different than what 
is expected from models without a density of states sin- 
gularity, where z = d or for those with straight Coulomb 
interation (z = 1 observed in the quantum Hall case). 
However, it is consistent with several field-theoretic sce- 
narios which give z = 2 (I116p . 




FIG. 48 The real part of universal scaling function E, as a 
function of a) = oj/T. There is a Drude-like peak from the 
inelastic scattering between thermally exciting carrier, which 
falls off at order w ~ 1. At larger lD there is a crossover to a 
coUisionless regime. Importantly the function gives different 
values in the lj and u ^ <x limits. 



I should mention that as presented above, frequency 
scaling does not necessarily give us any information on 
exponents etc. that we could not also get from the tem- 
perature scaling. Here the primary importance of fre- 
quency scaling was in its ability to demonstrate how char- 
acteristic time scales diverge at the QPT, which leaves 
the temperature itself as the only energy scale in the 
problem. However, as discussed by Damle and Sachdev 
(|117p . the fact that at the QPT itself, response func- 
tions can written as a universal function of lu/T, means 



that one does not necessarily expect the same behav- 
ior in the ui = 0, T ^ (incoherent) limit as in the 
uj ^ 0, T = (phase-coherent) limit . Since all DC 
experiments are in the former limit, while the vast ma- 
jority of theoretical predictions are in the latter, finite 
frequency measurements can in fact given unique insight. 
It was predicted by Damle and Sachdev that at the 2D 
superfluid-insulator transition the conductivity is equal 
to two different universal numbers of order e'^/h in the 
uj/T ^ and ui/T — > oo limits as shown in Fig. 48 This 



is true even if u> and T are both asymptotically small. 
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